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2 ■ ABSTRACT 

This paper revisits the classical Kennicutt method for inferring the stellar 
initial mass function (IMF) from the integrated light properties of galaxies. The 
• large size, uniform high quality data set from the Sloan Digital Sky Survey DR4 

O . is combined with more in depth modeling and quantitative statistical analysis to 

search for systematic IMF variations as a function of galaxy luminosity. Galaxy 
Ha equivalent widths are compared to a broadband color index to constrain 
the IMF. This parameter space is useful for breaking degeneracies which are 
O ■ traditionally problematic. Age and dust corrections are largely orthogonal to 

IMF variations. In addition the effects of metallicity and smooth star formation 
^ ' history e-folding times are small compared to IMF variations. We find that for 

■ the sample as a whole the best fitting IMF slope above 0.5 M0 is F = 1.4535 with 

a negligible random error of ±0.0004 and a systematic error of ±0.1. Galaxies 
brighter than around M^ 0.1 = ~20 (including galaxies like the Milky Way which 
has M^o.i ~ "21) are well fit by a universal F ~ 1.4 IMF, similar to Salpeter, 
and smooth, exponential star formation histories (SFH). Fainter galaxies prefer 
steeper IMFs and the quality of the fits reveal that for these galaxies a universal 
IMF with smooth SFHs is actually a poor assumption. Several sources of sample 
bias are ruled out as the cause of these luminosity dependent IMF variations. 
Analysis of bursting SFH models shows that an implausible coordination of burst 
times is required to fit a universal IMF to the M-rfi.i = —17 galaxies. This leads 
to the conclusions that the IMF in low luminosity galaxies has fewer massive 
stars, either by steeper slope or lower upper mass cutoff, and is not universal. 
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Subject headings: galaxies: evolution — galaxies: stellar content — stars: mass 
function 



Introduction 



A precise measurement of the stellar initial mass function (IMF) and its functional de- 
pendence on environmental conditions would impact astronomy over a wide range of physical 
scales. It would be of great help to theorists in untangling the mysteries of star formation 
and it is a key input in spectral synthesis models used to interpret the observed properties 
of galaxies both nearby and in the early universe. 

The current question is whether t he IMF is universal- the same regardless of time and 
environmental conditions. iKennicuttI (jl998al ) concisely states the current understanding of 
IMF universality. It is difficult to believe that the IMF is universal given the diversity 
of galaxy types, environments, star formation rates, and populations within galaxies over 
the range of observable lookback times. On the other hand, while IMF measurements do 
vary they are all consistent with a universal IMF within measurement errors and sampling 
statistics. The only way to proceed then is to strive for smaller measurement errors and 
improving sample sizes. 

A definitive theoretical derivation of the IMF does not yet exist. Theoretical approaches 
to the IMF usually center around the Jeans mass, Mj, the mass at which a homogeneous 
gas cloud becomes unstable. At first the collapse of a clo ud is isothe rmal a nd the Jeans 
mas s decreases which leads to fragmentation of the cloud (lHoylelll953l ). Both iReed (Il976l ) 
and ILow fc Lynden-Belll (Il976l ) suggested that at some point during the cloud collapse the 
line cooling opacity becomes high enough that the collapse is no longer isothermal. At this 
point the Jeans mass increases and fragmentation stops. The minimum Jeans mass is the 
smallest fragment size at this point and it provides a lower limit to size of the stars formed. 

Authors have calculated Jeans masses and minimum Jeans masse s using a variety of 
metho ds. In the classical derivation of the Jeans mass Mj oc T'^/^p"^/^. ILow fc Lynden-Bell 



( 19761 ) finds that Mj^min oc m~^^/^K~^/'' where m is the mass of gas atoms or molecules and 
K is the opacity at the final fragmentation. More recently turbulence in clouds has been 
studied. IPadoan. Nordlund. fc Joned (119971 ) found Mj^in oc n~^/'^T'^a~^ where n is the 
number density and cr^ is the velocity dispersion of the gas. Other investigators have looked 
at the hierarchical fractal geome try of mole c ular c louds, thought to aris e from turbulence, 
as a g enerator of the IMF (e.g. lElmegreenl (119971 )). On a related note lAdams &: Fatuzzo 
(119961 ) point out that molecular clouds exhibit structure on all resolvable spacial scales 
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suggesting that as no characteristic density exists for the clouds neither does a single Jeans 
mass. They develop a semi-empirical model for determining the final masses of stars from 
the initial conditions of molecular clouds without invoking Jeans mass arguments and use it 
to construct IMF models. The key components of their model are sound speed and rotation 
rate of cloud cores and the idea that stars help determine their final masses through winds 
and outflows. 

But from the beginning the study of the IMF has been driven by measurements. In 
1955 Salpeter was the first to make a meas urement of the IMF inferring it from his observed 
stellar luminosity function (jSalpeterlll955l ). We parameterize the IMF by: 

dn r -0.5 for 0.1 < m/Mo < 0.5 
dlogm °^ \ -r for 0.5 < m/M© < 120 ^ ' 

following iBaldry fc Glazebrookl (120031 ). Salpeter found that F = 1.35. It is often overlooked 
that his measurement only covered masses for which 0.4 < tti/Mq < 10. Nonetheless his 
original measurement is surprisingly consistent with modern values over a wide range of 
masses. The IMF in equation [1] is similar to the Salpeter IMF for F = 1.35. The difference 
is that there are fewer stars with masses less than 0.5 Mq. We adopt a two part power law 
as there is agreement amongst several authors that there is a change in the IMF slope near 
0.5 Mq (lKroupall200ll ). The technique we will use is not sensitive to the IMF at low masses 
so we assume a constant value in that regime. 

Salpeter's idea continues to be used today in IMF measurements of resolved stellar pop- 
ulations. The technique can be applied to field stars as well as clusters. However Salpeter's 
method has several inherent limitations. The nature of stars presents a challenge. On the 
high mass end stars are very luminous, but live only a few million years, while on the low 
end stars are faint but have lifetimes many times longer than the current age of the uni- 
verse. There are very few star clusters which are both young and close enough to allow us 
access to the IMF over the full mass range. In addition, the main sequence mass-luminosity 
relationship is a function of age, metallicity, and speed of rotation in addition to mass. It is 
not yet well-known at the low an d high mass extremes. Unresolved binaries can also affect 
the measured IMF (jKroupall200ll ). The light from unresolved binaries is dominated by the 
more massive of the pair. As a result the less massive star is typically not detected which 
l eads to a systematic un der-counting of low mass stars. 42% ± 9% of main sequence M s tars 
(iFischer fc Marcylll992l ) and 43% of main sequence G stars (jPuquennoy &: Mayor! Il99ll ) are 
primary stars in multiple star systems. These are both lower limits as some companions may 
have eluded detection. As roughly half of stars are in multiple systems it has a potentially 
large effect on the observed luminosity function. 



Except for at the high mass end field stars in the solar neighborhood offer the best 
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statistics for IMF measurements. However the solar neighborhood IMF is f ound to be de- 
ficient in r aassive stars when compared to other galaxies. For example the iMiller &: Scalo 
andlScald (I19861) solar neighb orhood IMFs are rejected by integrated lig ht approaches 



I.e. 



Kennicutt 



IlQsj), hereaf ter K83. lKennicutt. Tamblvn. fc Congdo3 Jl994h (KTC94), and 



Baldrv fc Glazebrookl (120031 ). 



Analysis of individual star clusters can be used to detect IMF variations. As methods 
and data quality can vary between aut hors comparisons between individual clusters are 
difficult. However iPhelps fc Janed (Il993l ) studied eight young open clusters with the same 
technique. On the extremes they measured F = 1.06ib0.05 for NGC 663 and F = 1.78 ±0.05 
for NGC 581 over masses from around 1 to 12 Mq. 



Kroupal (120011 ) considered the case in which one could have perfect knowledge of the 



masses of all stars in a cluster. Clusters have a finite size so even with no measurement errors 
uncertainty arises from sampling the underlying IMF. He shows that the observed scatter 
in IMF power law values above 1 Mq c an be accounted for b y sampling bias for clusters 
with between 10^ and 10^ members (the iPhelps &: JanesI (119931 ) clusters have memberships 
in this range). Furthermore, dynamical evolution of clusters can affect measurements of the 
IMF. This can happen by preferentially expelling lower mass sta r s from the cluster and by 
breaking up binary systems within the cluster. In total iKroupal (120011 ) finds that for stars 
with m > IM0 the spread in the observed IMF power-law slopes in clusters is around 1 
when both binary stars and sampling bias are considered even when the underlying IMFs 
are identical. 

Stochastic processes can also infiuence the ability to detect systematic IMF variations. O 
and B stars produce ionizing photons and cosmic rays which affect the surrounding nebula. 
However the probability of creating one of these massive stars is comparatively small. If 
one of these massive stars happens by chance to form first it may drive up the nebular 
temperature a nd depress the forma tion of less massive stars compared to regions without 
massive stars ( jRobberto et al.l 120041 ). 



An alternative to IMF measurements based on the stellar luminosity functions of re- 
solved stellar populations is to infer the characteristics of stellar populations from the inte- 
grated hght of galaxies using spectral synthesis models. The advantage of using integrated 
light techniques is that many of the problems plaguing IMF investigations of resolved stellar 
populations are avoided. Stochastic effects are washed out over a whole galaxy. Unresolved 
multiple star systems are irrelevant. The number of observable galaxies is large and their 
environments span a much larger range than those of Milky Way clusters. Integrated light 
techniques can be applied to the high redshift universe. This creates a strong motivation 
to develop IMF techniques and test them for galaxies in the local universe which can later 
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be used to probe earlier stages of galaxy formation. The assumption of a universal IMF 
has a huge influence on the interpretation of the high redshift universe. The reionization of 
the universe and the Madau plot- the global star formation rate of the universe over cos- 
mic history- are two areas where IMF variations could impact the current picture of galaxy 
formation and evolution. 

On the downside conclusions are dependent on the stellar evolution models used, which 
are not well constrained at high masses or with horizontal branch stars at low metallicities. 
The biggest problem is that changes in the IMF, star formation history (SFH) or age of a 
galaxy model can have similar effects in the resulting spectral energy distribution (SED). 
Any integrated light technique needs to address these degeneracies. It is also difficult to 
probe the IMF at sub-solar masses using integrated light. 

On the level of galaxies the concept of a universal IMF has recently become more com- 
plicated. A number of recent studies have investigated the effect of summing the IMF in 
individual clusters o ver a galaxy in the presence of power law star cluster mass functions. 



Weidner fc Kroupal (120051 ) argue that the integrated galaxial IMF will appear to vary as 



function of galactic stellar mass even if the stellar IMF is universal. However lElmegreen 



( 120061 ) argues that the galaxy wide IMF should not differ from the IMF in individual clus- 
ters based on analytical arguments and Monte Carlo simulations. Our approach can only 
measure the IMF averaged over whole galaxies and cannot address this distinction. Even so 
systematic variations of any kind from the Salpeter slope have not been measured outside of 
some evidence for non-standard IMFs in low su rface brightness galaxies (LSB) and galaxies 



experiencing powerful bursts of star formation (jElmegreeru |2006| ) . Either way, observational 



evidence for systematic variations of the IMF in galaxies would be highly valuable. 

The plan of this paper is as follows: §2 explains our method for constraining the IMF. 
In §3 we describe the SDSS data and our sample selection. §4 details our modehng scheme. 
In §5 we discuss our statistical techniques. §6 reports our results. In §7 we check our results 
against the B.Sa distribution of the data and §8 presents our conclusions. 



2. Methodology 

This paper revisits the "classic" method of K83 (and the subsequent extension KTC94) 
to constrain the IMF of integrated stellar populations. The method takes advantage of the 
sensitivity of the Ha equivalent width (EW) to the IMF. K83 showed that model IMF tracks 
can be differentiated in the {B — V) — log (Ha EW) plane. 

The total flux of a galaxy at 6565A is the combination of the underlying continuum flux 
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plus the flux contained in the Ha emission hne. The Ho; flux and the continuum flux have 
different physical origins, both of which can be used to gain physical insights into galaxies. 

In the absence of AGN activity the Ha flux is predominantly caused by massive O and 
B stars which emit ionizing photons in the ultraviolet. O and B stars are young and found in 
the regions of neutral hydrogen in which they formed. In Case B recombination, where it is 
assumed that these clouds are optically thick, any emitted Lyman photons are immediately 
reabsorbed. After several scattering events the Lyman photons are converted into lower 
series photons (including Ha) and two photon emission in the I'^S continuum. These 

photons experience smaller optical depths and can escape the cloud. The transition prob- 
abiUties are weakly dependent on electron density and temperature and can be calculated. 
Through this process the measured Ha flux can be converted into the number of O and B 
stars currently burning in an integrated stellar population. 

However Case B recombination is an idealized condition and it is possible that ionizing 
photons can escape the cloud without this processing, a situation known as Lyman leakage. 
As such the Ha flux is a lower limit on the number of O and B stars present. 

The continuum flux of a galaxy is due to the underlying stellar population. At 6565A 
the continuum is dominated by red giant stars in the 0.7-3 range while the Ha flux 
comes from stars more massive than 10 M©. 

The EW is defined as the width in angstroms of an imaginary box with a height equal 
to the continuum flux level surrounding an emission or absorption line which contains an 
area equal to the area contained in the line. This is effectively the ratio of the strength of a 
emission or absorption line to the strength of the continuum at the same wavelength. Given 
the physical origins of the Ha ffux and the continuum at 6565A the Ha EW is the ratio of 
massive O and B stars to stars around a solar mass. Therefore the Ha EW is sensitive to 
the IMF slope above around 1 M0 and can be used to probe the IMF in galaxies. 

As mentioned in the introduction several degeneracies plague the study of the IMF from 
the integrated hght properties of galaxies. Variations in the IMF, age, metallicity, and SFH 
of galaxy models can all yield similar effects in the resulting spectra. For example, increasing 
the fraction of massive stars, reducing the age of a galaxy, lowering the metallicity, and a 
recent increase in the star formation rate will all make a galaxy bluer. 

Metallicity effects were not discussed in either K83 or KTC94 and galaxy ages were 
assumed, 15 Gyr for K83 and 10 Gyr for KTC94. 

The SFH in K83 is addressed by calculating models with exponentially decreasing SFHs 
for a range of e-folding times, as well as a constant and a linearly increasing SFH. In the 
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[B — V) — log (Ha EW) plane the effect of varying the SFH e-folding time is orthogonal to 
IMF variations. However this is only true for smoothly varying exponential and linear SFHs. 
Discontinuities, either increases (bursts) or decreases (gasps), in the star formation rate can 
affect the Ha EW relative to the color in ways similar to a change in the IMF. Along with 
the exponential SFHs KTC94 also uses models with instantaneous bursts on top of constant 
SFHs. However this was done to access high EWs at an age of 10 Gyr rather than to fully 
flesh out the effects of SFH discontinuities. 

The assumption of smoothly varying SFHs is a key assumption in our analysis. Most 
late-type galaxies are thought to form stars at a fairly steady rate over much of r ecent time 



although bursts of star formation may play a significant role in low mass galaxies (IKennicutt 



1998bl ). For these galaxies smoothly varying SFHs are justified. However there are other 



galaxies clearly in the midst of a stron g burst of star formati on (e.g. M82) and dwarf galaxies 



with complex SFHs, e.g. NGC 1569 (lAngeretti et al.ll2005l ). for which this assumption is a 



poor one. The effects of violations of our SFH assumptions are described in detail in the 
results section. 



3. The Data 

K83 cites four major sources of error all of which are improved upon or eliminated by 
the high, uniform quality of SDSS spectroscopic and photometric data. 

The Ha fluxes in his sample can be contaminated by nonthermal nuclear emission. 
In the updated investigation, KTC94, this problem is addressed by removing objects with 
known Seyfert or LINER activity and luminous AGN. The SDSS spectra allow measure- 
ments of emission line ratios which ca n be used to separate star forming galaxies from AGN 



flBaldwin. Philhps. fc Terlevichlll98lr ) 



The second problem is that the Ha emission flux will be underestimated if the underlying 
st ellar absorption of Ha is n ot taken into consideration. The narrow band filter photometry 



of IKennicutt fc KentI (119831 ) could not measure this effect for individual galaxies so a fixed 



ratio was assumed for all galaxies. The SDSS spectroscopic pipeline does n ot take this into 



account either. We use the Ha fluxes measured from the SDSS spectra by iTremonti et al. 



(120041 ) which fit the continua with stellar population models to more accurately measure the 
Ha emission. While the SDSS pipeline method is sufficient for strong emission lines the Ha 
absorption EW can be as large as 5 A which is significant for weaker emission lines. 

Thirdly, their narrow band Ha imaging includes [N II] emission which is corrected for 
by assuming a constant [N II] /Ha ratio. The Ha and [N II] emission lines are resolved in 
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the SPSS spectra so the re is no need for a correction. This is a significant improvement. 



Kennicutt &: KentI (Il983l ) found from a hterature survey that the mean value of the Hq;/(Hq; 
+ [N II]) ratio is 0.75 ± 0.12 for spiral galaxies and 0.93 ± 0.05 in irregular galaxies. These 
mean corrections were applied uniformly to the K83 data. In KTC94 a uniform correction 
was applied using [N II] /Ha = 0.5. For comparison in our sample the mean value of the 
HQ;/(Ha + [N II]) ratio is a strikingly similar 0.752 and the mean [N Il]/Ha ratio is 0.340. 
However [N II] /Ha ranges from 0.0 to 0.6. If our Ha and [N II] lines were blended, applying 
a fixed correction would introduce errors of as much as 25% in the Ha EWs of individual 
galaxies. 

Lastly, in both K83 and KTC94 extinction corrections were addressed by plotting data 
alongside models which were either assumed an average value for the extinction or were 
unextincted. The Balmer decrement (Ha/H/3) can be measured from SDSS spectra which 
allows for extinction corrections for individual galaxies. 

Another major advantage of this study is that the sample size is much larger than that 
of K83 and KTC94. The KTC94 sample contains 210 galaxies, whereas ours has ~ 10^. 
The large sample size allows us to investigate IMF trends as functions of galaxy luminosity, 
redshift and aperture fraction with subsamples larger than the entire KTC94 sample. 

There is a key disadvantage to this method as well. K83 and KTC94 were able to adjust 
the sizes of their photometric apertures to contain the entire disk of individ ual galaxies to a 



limiting isophote of 25/i given by the RC2 catalog (iKennicutt fc Kentlll983l ). The advantage 
of narrow band measurements of Ha EW is that they can cover a much larger aperture and 
match the broadband measurements set to match the physical size of individual galaxies. 
The SDSS has fixed 3" spectroscopic apertures. This problem is partly offset by using 
matching 3" photometry apertures from the SDSS. However this introduces aperture effects 
as observed galaxies have a wide range of angular sizes due to the range of physical sizes 
and distances present in the local un iverse. This is significant as radial metallicity gradients 



[e.g. IVila-Costas fc Edmunds! (119921 )) are observed in spiral galaxies which could incorrectly 
be interpreted as radial IMF gradients. 

Even so, this method can constrain the IMF within the SDSS apertures. For our program 
galaxies 23% of the total light falls in the SDSS apertures. The fact that the more distant 
galaxies are more luminous and tend to be larger helps to balance out the larger physical 
scales of the fixed aperture size at greater distances. On average 17% of the light falls in 
the aperture for the faintest galaxies while it is 25% for the brightest bin. In spite of their 
limited size the SDSS apertures still contain a great diversity of stellar populations which 
make this data set an excellent test bed for IMF universality. We will present extensive tests 
of aperture effects below. 
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3.1. Sample Selection 



The sample is selected from Sloan Digital Sky Survey data. The project goal of the 
SDSS is to image a qu arter of the sky in five optical bands with a dedicated 2.5 m tele- 
scope ( lYork et al.l 120001 ). From the imaging 10^ galaxies and 10^ quasars will be selected 
for spectroscop i c foU owup. Photometry is done in the ugriz fi lter system described by 
Fukugita et al.l (Il996l ). Magnitudes are on the arcsinh system ( iLupton. Gunn. fc Szalay 
I999I ) which approaches the AB system with increasing brightness. Spectra are taken 
with a mul ti-object fiber spect rograph with wavelength coverage from 3800A to 9200A and 
R ~ 1800 JUomoto et al.lll999h . 



Our sample is a sub-sample of the Main Galaxy Sample from SDSS DR4 (lAdelman-McCarthy et al 



20061 ). The Main Galaxy Sam ple (MGS) targets galaxies with r < 17.77 in Petrosian magni- 
tudes ( IStoughton et al. I l2002h . All galaxies in the MGS are strong detections so the differences 
between luptitudes and the AB system can be ignored. In order to avoid fiber crosstalk in 
the camera an upper brightness limit of g = 15.0, r = 15.0 and i = 14.5 is imposed. Targets 
are selected as galaxies from the imaging by comparing their PSF magnitudes to their de 
Vaccouleur's and exponential profile magnitudes. Exposure times for spectroscopy are set 
so that the cumulative median signal-to-noise satisfies (S/N)^ > 15 at (7 = 20.2 and i = 19.9 
fiber magnitudes. The time to achieve this depends on observing conditions but always 
involves at minimum three 15 minute exposures. 

Due to the construction of the spectrograph fibers cannot be placed closer than 55" to 
each other. This may be a source of bias in the sample. Cluster galaxies may be preferentially 
excluded fror n the sample. The SDSS collabo ration has plans to quantify this effect in the 
near future (lAdelman-McCarthy et al.l l2006l ). LSBs are excluded from the MGS with a 



surface brightness cut which may also bias the sample (IStoughton et al.l 120021 ) . There is 



some evidence that LSBs may have IMFs which differ from a universal IMF. iLee et al. 



( I2OO4J ) find that the comparatively high mass-to-light ratios of LSBs can be explained with 
an IMF deficient in massive stars relative to normal galaxies. 



^Our sample begins with the fourth data release (DR4) of the SDSS (lAdelman-McCarthy et al. 

20061 ). DR4 covers 4783 deg^ in spectroscopy for a total of 673,280 spectra, 567,486 of 



which are galaxy spectra. 429,748 of these have fiags set indicating they are part of the 
MGS. DR4 also includes special spectroscopic observations of the Southern Stripe which 
were not selected by the standard algorithm but which nonetheless have the TARGET_GALAXY 
fiag set in primTarget which usually indicates membership in the MGS. These objects are 
i dentified by comparing their s pectroscopic plate number to the list of special plates in 
( lAdelman-McCarthy et al.ll2006l ) and rejected. This leaves 423,285 spectra. 
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The first round of cuts to our sample address general data quality. While the overall 
quality of SDSS data is high there are a handful of objects with pathological values for 
one or more parameters. First we require that all parameters of interest have reasonable, 
real values. This means that the Petrosian and fiber magnitudes must be between and 
25 and have errors smaller than 0.5 in all five ugriz bands. Line fiux errors are capped at 
10~^^ erg cm~^ s~^ and equivalent widths at lO^A. For most of the parameters less than 
1% of objects fail this loose requirement. However 6.6% of the objects fail the Petrosian 
magnitude error requirement. This is most likely due to the photometric pipeline having 
a difficult time defining the Petrosian radius. As such this constraint is potentially biased 
against LSBs or galaxies with unusual morphologies. In defense of this cut we later bin 
our data by luminosity and aperture fraction both of which are determined in part by the 
Petrosian magnitudes and also by i^'-corrections determined from them. In addition, limiting 
flux errors to 50% is hardly unreasonable. Altogether 391,160 galaxies pass these combined 
requirements, which is 92.4% of the MGS. 

Galaxies from photometry Run 1659 are removed because of a known problem with the 
photometry. This excludes 4,485 galaxies (1.1%) from a continuous strip on the sky and 
should not be a source of bias. We place a further constraint on the z band fiber magnitude 
requiring that the error be less than 0.15. The z band generally suffers from the most noise 
so this requirement ensures that the fiber magnitude quality is good enough to minimize the 
chance of erroneous if-corrections which could affect our colors. Only 2,866 (0.7%) MGS 
galaxies fail this test. 

Combining the general data quality requirements, the Run 1659 rejection and the fiber 
z band error limit leaves 386,647 galaxies (91.3% of the MGS). 

The next round of cuts to our sample, while necessary, have clear astrophysical impli- 
cations. Many of the objects in the MGS have AGN components. As we are interested in 
studying only the underlying st ellar populations of these objec t s AG N must be removed. 



This is done using the classical [Baldwin. Phillips. Sz Terlevichl ( 119811 ) diagram comparing 



the logarithms of the [O III A5007]/H/j and [N II A6584]/Ha emission line ratios. We used 
the criterion of Kauffmann et al. J2003h where objects for which 



■°g<l°"'l/"'^)> los(lNII]/Hl)-0.06 + ^-^ <2) 



are classified as AGN and rejected. Following iBrinchmann et al.l (12004J ) we require the S/N 



of the Ha, H/5, [O III] and [N II] lines to be at least 3 to properly classify a galaxy as a 
star forming one. 131,807 galaxies (34.1% of MGS objects surviving our first round of cuts) 
survive this cut. 



The above cut automatically rejects any galaxies with weak [O III] and [N II] lines. This 
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excludes galaxies with weak star formation. To a lesser extent metal poor g alaxies which also 



have weak [N II] emission are rejected as well. iBrinchmann et al.l (120041 ) also define a low 
S/N star forming class of galaxies which we identify and keep in our sample. These are the 
galaxies which have not already been classified as star forming or AGN by strong lines and 
equation [21 nor have been identified as low S/N AGN by [N II A6584]/Ha > 0.6 with S/N 
> 3 in both lines, yet still have Ha with S/N at least 2. 79,548 (20.6%) of the sample falls 
into the low S/N star forming galaxy category. Combining the two classes 211,355 (54.7%) 
of the galaxies survive the AGN cut. 

The AGN cut also has a strong luminosity bias for two reasons. Galaxies with AGN 
compo nents tend to be b righter. The bimodal distribution of galaxies in color-magnitude 



space (IBaldry et al.ll2004l ) also plays a role. The most luminous galaxies are predominantly 
red with minimal star formation and thus weak emission lines. Luminous galaxies are rejected 
both for having AGN components and for having low S/N emission lines. Over 95% of 
galaxies fainter than Mr,o.i = —19 meet this criteria, but by M,.^o.i = —24 the fraction is only 
38.9%. 

A color bias is also introduced by the AGN cut. Over 95% of galaxies bluer than 
{g — r)o.i = 0.6 pass, which drops to 24% by {g — r)o.i = 1.2. This is mainly due to the S/N 
requirement for the emission lines. Redder galaxies tend to have weak emission lines and are 
rejected. 

The Balmer decrement is used for the extinction correction so Ha and H/5 S/N are 
required to be at least 5 to reduce errors. 214,912 galaxies (55.6%) have H/? S/N > 5 which 
is the more restrictive of the two criteria. This cut has a clear luminosity bias. Roughly 
85% of galaxies with Mr,o.i > —20 satisfy this requirement, but this fraction decreases with 
increasing luminosity until only 12.9% survive at M^o.i = —24. This is again due to the 
bimodal distribution of galaxies. The luminous red galaxies with weak emission lines are 
rejected. 

There is also a color bias. Over 93% of the bluest galaxies blueward of {g — r)o.i = 0.8 
survive the cut while only 23% of the reddest pass this requirement. As previously mentioned, 
by nature the reddest galaxies have weak Balmer lines as a result of their low SFRs and are 
preferentially rejected. 

A redshift cut of 0.005 < z < 0.25 is applied to ensure that peculiar velocities do not 
dominate at low redshift and to limit the range of galaxy ages. 384,349 galaxies (99.4%) 
meet this criteria. Nearly all galaxies from M^^o.i = —17 to —23 survive this cut. On the 
low luminosity end only 27.4% of M^^o.i = —14 galaxies are distant enough to pass and on 
the high end 93.1% of M^^o.i = —24 galaxies are close enough to survive. Only 63% of the 
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bluest galaxies pass. Many of the blue galaxies that fail are actually H II regions of Local 
Group galaxies which are treated as their own objects by the SDSS pipeline so removing 
them actually improves the integrity of our sample. 

The stellar populations of galactic bulges can be significantly different from those in 
the spiral arms. The SDSS fibers have a fixed aperture of 3" so over the large range of 
luminosities and distances in the MGS aperture affects can become very important. To 
remove outliers we require that at least 10% of the light from a galaxy falls within the 
spectroscopic aperture. This is done by comparing the Petrosian magnitude to a fixed 3" 
aperture fiber magnitude. Both of these quantities are calculated for all objects in the 
SDSS by the photometric pipeline. 371,777 (96.2%) galaxies survive this cut. This cut 
rejects proportionally more faint galaxies; 99.2% pass at M^^o.! = ~24 compared to 50.9% at 
Mr-,0.1 = —14. The aperture cut has low sensitivity to color. 

The intersection of the AGN, Balmer line S/N, redshift and aperture fraction cuts leaves 
140,598 galaxies- 36.4% of the high quality MGS data defined by our first round of cuts and 
33.2% of the MGS as a whole. 

At this point three final cuts are applied to the sample. One galaxy is removed for 
surviving all criteria, but having a negative Ha EW. 

The extinction of individual galaxies is estimated using the Balmer decrement (Ha/H/5 
emission flux ratio) for each galaxy. Given Case B recombination, a gas tem perature of 
10,00 0K and density of 100 cm~^ the Balmer decrement is predicted to be 2.86 ( Osterbrock 



989f) . This ratio is weakly dependent on nebular temperature and density. lOsterbrock 
(119891 ) lists values down to 2.74 for Case B recombination in environments where both the 
temperature and electron density are high. 3.2% of the galaxies suffer from the problem 
that the Balmer decrement is less than 2.86 and 2.1% have a Balmer decrement below 2.74. 
538 galaxies (0.4%) have Balmer decrements more than 3(j below 2.74, which is around 6 
times more than expected. This is does not suggest a problem with the Case B assumption 
as the predicted Balmer decrements for Case A recombination are nearly identical in each 
temperature regime. 

To understand the reason behind this problem around 100 of the offending spectra were 
inspected revealing a few different causes for the problem. Around 80 galaxies are at redshifts 
where the telluric O I A5577 line affects the measurement of H/3. Many of these galaxies have 
very strong emissio n lines with extr e mely weak stellar absorption. Using the SDSS pipeline 



values instead of the lTremonti et al.l (120041 ) values yields acceptable Balmer decrements. The 
rest are galaxies with low flux where the Balmer lines are in absorption. This shows there 
are a few cases where attempting to fit the underlying stellar absorption lines fails and this 
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failure is not reflected in the error values. These galaxies are rejected without any apparent 
introduction of bias. 

All of these cuts combined leave 140,060 galaxies, which is 33% of the MGS and 36% of 
the high quality MGS data. After removing duplicate observations there are 130,602 galaxies 
in our sample. Of these objects 1.7% overlap with the Luminous Red Galaxy Sample. 

Overall the bulk of the galaxies are removed by AGN rejection and the H/5 S/N require- 
ment, with the rest of the cuts having little effect. Both of these cuts are necessary. AGN 
must be removed to ensure that the Ha emission represents the underlying stellar population 
and not an accretion disk. Our method requires that the galaxies have measurable Balmer 
emission lines. This coupled with the need for accurate extinction corrections justifies the 
Balmer line S/N requirement. 

Our cuts bias our sample by preferentially excluding galaxies at both luminosity and 
both color extremes. The faintest galaxies are most affected by the redshift and aperture cuts 
while the luminous galaxies succumb to the H/3 S/N requirement. At the red extremes it is 
the H/3 S/N and AGN requirements that play equally large roles, while the bluest objects are 
primarily rejected by the Hubble flow redshift requirement. Although our sample is biased 
by our cuts each one is a necessary evil. We do not attempt to correct this bias, but we 
remind the reader that the following results are only representative of actively star forming 
galaxies without any AGN activity. 

The aim of this paper, however, is to test IMF universality. If the IMF is truly universal 
it should be universal in any subsample of galaxies. The fact that our sample is slightly biased 
with respect to luminosity and color is not a significant barrier to achieving our goal. 



3.2. Corrections 



The SDSS includes a number of different calculated magnitudes. We use the fiber mag- 
nitudes which are 3" fixed aperture magnitudes. Although it was not the c ase in earlier ver 



sions of the photometric pipeline, fiber magnitudes are now seeing corrected (lAbazajian et al. 



20041 ). The fiber magnitudes were not originally intended for science purposes but rather to 
get an idea of how bright an object will appear in the spectrograph. We use the fiber mag- 
nitudes to reduce the aperture effects arising from comparing a 3" spectroscopic aperture to 
Petrosian magnitudes. Originally the SDSS used "smear" exposures to correct spectra for 
l ight falling outside th e 3" aperture due seeing, guiding errors and atmospheric refraction 
(jStoughton et al.ll2002l ). The smear technique was later found to be an improv ement only for 
high S/N point sources and its use was discontinued (lAbazajian et al.ll2004j ). The spectra 
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here are not seeing corrected. 

After paring the sample to its final size a number of corrections must be made to both 
the photometric and spectroscopic data. G alactic reddening from the Mi l ky W ay must be 
corrected for. SDSS database includes the ISchlegel. Finkbeiner. fc Davis (|l998[ l dust map 
values for each photometry object. 

The extinction of individual galaxies is estimated using the Balmer decrement for each 
galaxy. The data is corrected assumin g th at 2.86 is the true value of the Balmer decrement 
using the Milky Way dust models of iPei Jl992h . The assumption of Milky Way dust is 



not significant as models of the dust attenuation in the Milky Way, SMC and LMC are 
nearly identical in the g band and redward. As aforementioned a few percent of our galaxies 
have Balmer decrements below 2.86. Our solution is to set the emission line extinction to 
A-v,i = 0.01 magnitudes for these galaxies. 

Massive young stars and their surrounding ionized nebulae tend to be embedded in their 
star forming regions more so than older, lower mass stars which have had time to migrate 
from their birth regions. As such nebular emission lines will ex perience more extinction than 



the continuum. ICalzetti. Kinney, fc Storchi-BergmannI (119941 ) find the ratio of emission to 



continuum line extinction is / = 2.0 ± 0.4. We assume this value to be 2.0 and correct the 
continuum and emission lines separately. This is the same value used by K83 and KTC94 
in their extinction corrected models. We note that the spatial geometry of the dust can 
influence the extinction law, but this complication is beyond the scope of this paper. 

Galaxy photome try is i^'-corrected to z = 0.10 using version 4.1.4 of the code of 



Blanton et al.l (j2003al ). This redshift is roughly the median of the sample and is selected 
to minimize errors introduced by the fC-corrections. The {g — r)o.i colors we use are the 
g — r colors we would observe if the galaxies were all located at 2; = 0.1. 

Stated explicitly the {g — r)o.i color is 

c={g- r)o.i = {g-kg-kg- 1.153Av,z//) - {r - h, - K - 0.834AvV/) (3) 

where kg and are ii"-corrections, A^ and A,, are Milky Way reddening values, and 1.153 
and 0.834 relate the V band extinction to the g and r bands assuming a Milky Way dust 
model. The corrected equivalent width is obtained as follows 



w 



Wo [(1 + Z)X 10-0-4(0.775Av,){l-l//)] -1 (4) 



where wq is the measured, uncorrected equivalent width. The 1 + z arises from the fact that 
the total flux in the Ha line is not affected by cosmological expansion of the universe but the 
flux per unit wavelength of the continuum is depressed by a factor of 1 + z. The following 
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term is the extinction correction. The 0.775 relates the V band extinction to the Ha hne 
assuming a Milky Way dust model and the 1 — 1// is due to the fact that the emission line 
and continuum experience different amounts of extinction as previously explained. 

Figure [1] shows the distribution of the galaxies in the color vs. Ha EW plane. 



3.3. Errors 

In order to conduct a likelihood analysis we need error estimates which take into account 
both the errors induced by the photometry and spectroscopy and those by the aforementioned 
corrections. The error in the corrected color, ctc, is given by 



cTc = 0.03 + 




+ (0.0440A,)' + ol (5) 



where cXg and are the Poisson errors in the observed g and r band photometry, / is the 
ratio of the emission line to continuum extinction, and is the error introduced by the K- 
corrections. The terms inside the square root in equation [S] are obtained through a straight 
error propagation of equation [3l The 0.03 outside t he square root is due to the systematic 



zero p oint errors of the SDSS filter system. Following lCalzetti. Kinney, fc Storchi-Bergmann 



( I1994J ) / is fixed at 2.0 and aj is set to 0.4. The error in emission line Ay is given by 
which is dependent on the fr a ction al uncertainty of the Ha and H/3 emission line fluxes. 



Schlegel. Finkbeiner. fc David (119981 ) reports a 16% error in their Milky Way reddening val- 
ues. Because a dust model is assumed reddening in the g and r bands are linearly related, so 
the error in — is a function of A^. This relationship combined with the 16% error yields 
the 0.0440 in equation [51 The median value of g band reddening is 0.10 so the errors intro- 
duced by the MW reddening correction are insignificant for the majority of objects. Errors 
in the redshift determination are negligible, typically 0.01%. The value of cr^ is estimated to 
be 0.02 by visual inspection of a plot of kg — fc^ as a function of redshift. 

For a typical galaxy the term involving cxj is the largest contributor to the extinction 
corrected color error. The median Poisson error from the photometry is 0.01 in both bands. 
The median value of cXc is 0.085. 
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The error in the corrected Ha equivalent width, a^, is given by 

2 2 

a^ = w 

Equation [7] is the result of propagating the errors in equation HJ neglecting the insignificant 
redshift errors. Again, the term involving o"/ is the largest contributor to the error for typical 
galaxies. The median error in the equivalent width is 17%. The median error bars for the 
sample are shown in Figure [H For comparison, K83 reports equivalent width errors of around 
10%, but the extinction is uncertain at the 20-30% level. 



^ 1 +0.5095 

Wo 




OfA 



MI 



(7) 



4. Models 

Model gal axy spectra were calculated usin g the publicly available PEGASE.2 spectral 



synthesis code (IFioc fc Rocca-Volmerangdll997l ). Models are calculated for ages from 1 Myr 



to 13 Gyr. 25 smoothly varying SFHs generated by analytic formulae are considered. The 
SFHs range from 19 exponentially decaying SFHs with time constants from 1.1 to 35 Gyr, a 
constant SFR, and four increasing SFHs which are proportional to 1 — exp^*/"^ where r is the 
time constant. The precise values of the time constants were selected to smoothly sample 
the Ha EW vs. {g — r)o.i plane. The metallicity of the stars is assumed to be constant 
with respect to time and calculated for Z = 0.005, 0.010, 0.020 and 0.025. Galactic winds, 
galactic infall and dust extinction are turned off. The dust extinction is not modeled because 
we have applied an extinction correction to the data. Nebular emission is calculated from 
the strength of the Lyman continuum. Emission line ratios are fixed. The model spectra are 
redshifted to 2; = 0.1 to match the redshift range of the data. 

The model parameter of interest is the IMF. Implicit in equation [T] is the assumption 
that the IMF does not vary as a function of time. 

The continua of galaxies are weakly infiuenced by low mass stars in the optical. This 
method is sensitive to the IMF for masses above around 1 Mq so the slope is fixed below 0.5 
Mq. Above 0.5 M0 models are calculated for 1.00 < F < 2.00, where F is incremented by 
0.05 between models. 

We treat our IMF model as though it has only one degree of freedom- F above 0.5 Mq. 
In truth it has three more as written: the lower and upper mass cutoffs and the point at 
which the slope changes. So how well are our assumptions justified? 



We have parameterized the IMF as a piecewise power law with t wo comp o nents . Piece- 
wise power laws are motivated by empirical fits to data starting with ISalpeterl (jl955l ). which 
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had only one component. By contrast the power law formulation of the IScald (119861 ) IMF has 
24 components. The log-normal distribut ion is a more physical choice as it can arise from 
stochastic processes. iMiller fc Scald (Il979l ) were the first to fit an observational measurement 
of the IMF with a log-normal distribution. The log-normal distribution is normalizable as it 
goes to zero smoothly at both extremes without any awkward truncation. Its main drawback 
is that it cannot fit any structure in the IMF over small mass ranges. 

Log-normal distributions have three degrees of freedom. This is less than the four our 
model has. However, we are not sensitive to IMF over the full range of masses which makes 
it much more difficult to fit the parameters of the log-normal distribution. Instead we use 
this piecewise model and lower the degrees of freedom through physical arguments. Many 
investigators find a change in slope in the IMF around 0.5 Mq. Our fixed lower end of the 
IMF is designed to be consistent to this. As our technique is not sensitive to this regime this 
assumption does not impact the results. 

The IMF must be normalizable because the total mass of a stellar population is finite. In 
our parameterization this is achieved by truncation at 0.1 and 120 Mq. This seems unphysical 
as the existence of brown dwarfs suggests that the IMF should continue below the hydrogen 
burning limit. However, stars at 0.1 Mq do not contribute much to the integrated light of 
galaxies. As we are not sensitive to stars in this mass range this choice is not unreasonable. 
In fact, truncating the IMF at 0.5 M0 yields models which are at worse differ by 0.002 in 
{9 ~ and 3% in Ha EW from those with low mass stars. The slope below 0.5 Mq has 
essentially no effect on our results. Only when the IMF is truncated at 0.9 Mq do the models 
differ at the level of the errors in the data. 

On the high mass end the choice of limit does matter. There is a physical upper limit 
to the size of stars associated with the Eddington limit. The value of this theoretical limit 
is not widely agreed upon. The largest stellar rn as s measured reliably , via a nalysis of a 
binary system, is 83 ± 5 Mq (IBonanos et al.l 120041 ). IWeidner &: Kroupal (120041 ) argue that 
given the large mass and youth of the star forming cluster R136 in the Large Magellanic 
Cloud stars in excess of 750 Mq should be present given a Salpeter IMF with no upper 
mass limit to stars, whereas no stars above 150 Mq are observed. An analysis of the Arches 
Cluster, the youngest observable cluster, gives an upper limit of 150 Mq based on Monte 
Carlo simulations although stars above 130 Mq are not detected (lFigeiil2005l ). The PEGASE 
model tracks only extend up to 120 Mq so this is the cutoff used. Another issue is that the 
physics and evolution of such high mass stars is not well known so the models themselves 
may be a significant source of error in this regime. 

The left half of Figure [2] shows the effects of varying the high mass cutoff in the IMF. 
The effect in the {g — r)o.i-Ha; EW plane is seen to be very similar to increasing the value 
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of r. Lowering the upper mass cutoff from 120 to 90 Mq has roughly the same effect as 
increasing F from 1.35 to 1.45. The relationship between the change in the upper mass cutoff 
(from 120 Mq) and the apparent change in F is AF ~ 0.005AMup and is roughly linear for 
upper mass cutoffs down to 50 Mq. The coefficient in the relationship is a week function 
of the age of the population ranging from 0.004 for 13 Gyr old populations to 0.006 for 100 
Myr old populations. 

The right half of Figure [2] shows the affect of adding a second break in the IMF at 10 Mq. 
Reducing the value of F over the 0.5-10 Mq range while keeping it fixed at F = 1.35 above 
10 Mq has a similar effect to decreasing F in a two component model. This illustrates one of 
the limitations of this model. At this point it is not possible to detect fine structure in the 
IMF slope or to state precise values for the IMF slope. In this limited space of observables 
the IMF models themselves are degenerate. What it does provide is a framework with which 
to detect variations in the IMF. Although we can construct similar tracks from different IMF 
models, we can still detect the differences between two groups of galaxies. 

While we will report our results as a function of F it must always be kept in mind that it 
is degenerate with the upper mass cutoff and other fine structure in the IMF at high stellar 
masses. 

The assumption of smoothly varying SFHs is of great consequence. In the event that 
a galaxy is experiencing or has recently experienced a burst our SFH assumption can lead 
to measured F values that are off by as much as 0.5. The effects of bursts are more closely 
examined in a later section. 

Within the assumption of smoothly varying SFHs much can be said about the effects 
of the IMF, metallicity, age, and SFH in the color-Ha EW plane. Figure [3] demonstrates 
these relationships. In both panels the ages of the models decrease from the upper left to 
lower right. The effects of the age of a stellar population are largely orthogonal to those of 
IMF variations. In Figure [3^ the effects of changing the functional form of the smoothly 
varying SFH with fixed metallicity and F are shown. SFH variation is degenerate with the 
IMF. However the effect is relatively small over a wide range of SFHs. The solid lines are 
exponentially decreasing SFHs with r = 1.1 Gyr where the bulk of the star formation occurs 
early in the galaxy's life. The dashed lines have SFHs that are increasing with time where 
most star formation occurs cit Icitc ages. The effect of variations in the form of smooth 
SFHs is larger at later ages and higher values of F but does not dominate the effects of 
IMF variations. With all other parameters fixed the range of smooth SFHs cause systematic 
uncertainties at the level of ±0.1 in F. 

In Figure [3j3 the effects of metallicity variations with fixed SFH and F are shown. The 
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metallicity variations are also degenerate with IMF variations. With all other parameters 
fixed between colors of 0.1 < {g — r)o.i < 0.4 the systematic uncertainty due to metallicity 
is less than 0.05 in F. This uncertainty increases to 0.35 at {g — r)o.i = —0.2 and 0.7. 

As aforementioned, the extinction correction is another potential problem. The arrows 
in Figure [3] show the length and direction of the extinction correction for typical galaxies in 
our sample. It is assumed that / = 2, but / = 1 and / = 4 are also plotted to show the 
potential effect of variations in /. The reddening vectors for / = 2 and 4 are fortuitously 
orthogonal to the IMF variations. Only when the continuum and emission extinctions are 
equal, when / = 1, do the extinction correction and variations in / become a larger con- 
cern than metallicity and SFH. However suc h low / ratios are not observed in galaxies 



(jCalzetti. Kinney, fc Storchi-Bergmanru (119941 ) and section 16.2.71) 



Figure m shows all 18,480 model points with F = 1.35. Models are interpolated in SFH 
history for fuller coverage of the color-Ha EW plane. For each value of F the models cover 
a stripe rather than a single line. It can be seen in the lower right of figure H] that the model 
become degenerate in F for old, red galaxies with weak current star formation. 



5. Statistical Techniques 

The data and models are compared using a "pseudo-x^" minimization. For various 
reasons (detailed below) the classical estimator assumptions are violated so we can not 
use traditional tables for error estimates but we can still use the as a statistical estimator 
as long as the confidence regions are calibrated by Monte Carlo (MC) techniques as we will 
do. We proceed as follows: for each galaxy i we have a measured {g — r)o.i color, Cj, and an 
Ha EW, Wi, and measurement errors and a^., given by equations [5] and [3 We also have 
model values c(F, Z, r, ip) and w(T, Z, r, ip) for a range of IMF slopes F, metallicities Z, ages 
t and SFHs ip. We can then construct a value as 

which is calculated by brute force. The goal of this paper however is to investigate the 
IMF with relatively simple measurements of the Ha EW and a broadband color. While 
making crude measurements of the mean stellar metallicities of individual galaxies is possible, 
disentangling age and SFH effects on an individual basis is a daunting task. Assuming that 
it is possible to do, it does not scale up well to the high redshift universe where observations 
will be of lower quality. 

It does not make sense to minimize x^ over all galaxies for a particular set of (F, Z, t, ip) 
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because we have no a priori reason to think that all of the galaxies should have the same 
metallicity or SFH. In fact we expect that they would not. The solution is to marginalize 
over metallicity, age and SFH for each galaxy such that 

x^(r) = min[x?(r,z,t,v^)]^_^^^ (9) 

This is somewhat unorthodox, because for some galaxies the data points are over-fitted, 
i.e. there will be a stripe in {c,w) space corresponding to a given T and we can get XiiX) 
values very close to zero (but not exactly because of the discrete nature of the model grid) for 
galaxies within the stripe. We note we also have partial degeneracy between parameters such 
as age and metallicity — they both shift the tracks in similar directions largely orthogonal 
to r (though not completely which is why we have a stripe in parameter space not a line). 
This makes it difficult to calculate the traditional "number of degrees of freedom." Despite 
these limitations it is clear that galaxies inconsistent with a particular F will still have large 
values of xli^)^ fo^' example a very blue galaxy with a low equivalent width in Figure |H 
The complication is that the stripes for similar F values overlap, and for red galaxies with 
low equivalent widths the stripes for vastly different F values overlap. As such, the IMF 
for an individual galaxy is only broadly constrained. Measuring a precise best IMF for 
an individual galaxy boils down to random chance and the discrete nature of the models. 
However, by summing xf over many galaxies the IMF is narrowly constrained for the sample 
being summed over as long as we are careful in our confidence region analysis. 

Because of this over-fitting and partial degeneracy we can not apply the textbook notions 
of the distribution, calculate degrees of freedom and choose Ax^ contours for different 
confidence regions. Further to this q and Wi are not truly independent variables. Both the 
colors and EWs are subject to the same extinction and reddening corrections which tie the 
errors together. For galaxies with z < 0.04 the Ha line is in the observed r band, although 
this only affects a relatively small number of galaxies in the sample almost all of which have 
MrfiA > —20. Also the direct statistical interpretation of is predicated on the assumption 
of normal errors. Equations O [6] and [7] reveal that our errors are complicated mixtures of 
individual measurements which are most likely Poisson distributed. Thus cTc and are 
unlikely to be normally distributed. Bursty SFHs can potentially create outliers which are 
statistically significant due to the fact that neither ctc or include a term for this difficult 
to quantify effect. The problem is even worse if the errors are non-symmetric which could 
potentially arise from the aforementioned bursty SFHs. In the case of non-symmetric errors 
the best value of F could erroneously be pulled away from the true value. 

Given all this we abandon the direct statistical interpretation of x^ and regard it as an 
estimator of the goodness of fit whose confidence re gions have to be ca librated empirically. 
We do this via MC simulations (as recommended by Press et al. ( 19921 )) where we simulate 
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data points for a given T with the correct error distributions and propagate everything 
through the analysis in the same way as for the actual data. For each of our MC simulations 
we add Poisson errors to the ugriz fiber magnitudes, Ha and H/3 fluxes and the observed Ha 
EW. We assume that these observed quantities have Poisson dominated errors- as members of 
the MGS they are high S/N measurements. The entire analysis described above is repeated, 
including a new extinction measurement and a recalculation of the i^-corrections. For each 
value of r we run 100 MC simulations to estimate the 95% confidence interval. Setting up the 
MC architecture in this way has the further advantage that we can use the same machinery 
to test the effect of systematic errors such as the violation of our smooth SFH assumptions, 
as we will do later. The main downside of course is that this approach is computationally 
intense. Run times for the 100 MC simulations are typically 18 days on a 2 Ghz desktop PC 
for the samples considered here. 

In practice it turns out that x^(r) is still a smooth well-behaved function with, not 
surprisingly, a quadratic minimum which has the advantage that we can then interpolate it 
to increase the resolution in the best-fitting F without incurring the additional computational 
expense. This arises of course from the fact that our estimator is similar to a traditional 
and is a good reason to stick with this similarity over some more exotic goodness of fit 
measure. An estimate of the systematic errors is discussed later. 

Regardless of how poorly a sample is modeled by a universal IMF the above method 
will still find a best fitting F and corresponding confidence region. We still expect to be 
small for a model that is a good fit and large for one that is not. One nuance in comparing 
between different sub-samples, as we will do, is that the samples are often of considerably 
different sizes. Because of this we choose instead to use the mean x^, instead, as a sample 
metric. This has the advantage that absolute values and confidence regions are more 
similar between the sub-samples, though we note that the confidence regions on x^ are still 
determined directly from our MC simulations. 



6. Monte Carlo Results 

Figure [5] shows the results of our analysis for the full sample of galaxies using just the 
observed data set. The "X" marks the best fitting IMF, where F = 1.4411 and x^ = 60549.7 
with x^ = 0.46. At F = 1.00 x^ = 3.61 while steeper IMFs are more heavily rejected with 
p = 9.00 at F = 2.00. 



For comparison several "classic" IMFs are also plotted in figure [5] at their approxi- 
mate equivalent values of F. With x^ = 7.73 the Miller fc Scalol ( 1979 ) solar neighborhood 
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IMF is a particularly bad fit. Two more recent solar neighborhood IMFs, IScald (119861 ) 
and iKroupa. Tout. &: Gilmord (Il993l ) , yield = 2.38 and 1.98 respectively . These results 
reinforce the conclusions of K83, KTC94 and iBaldry &: Glazebrookl (120031 ) that the solar 
neighborhood is not representative of galaxies on the whole as far as the IMF is concerned. 



On the other hand the I Scald (119981 ) IMF, established from a review of star cluster IMF 
studies in the literature, is a better fit than the best F value in our parameterization with 
= 0.43. This result highlights the degeneracy of the IMF models themselves in the color- 
Ha EW plane- two considerably different IMFs (one with one break and the other with two) 
fit nearly equally as well. 

The results of the MC simulation show that the 95% confidence region is 1.4432 < F < 
1.4443 for the data set as a whole. The MC simulation shows that additional data will 
not improve the overall results as the random errors are already small. Clearly and not 
surprisingly systematic errors, which are discussed later, dominate. 

The overall result of F = 1.4437 ± 0.0005 is steeper tha r i the original Salpeter value of 
F = 1.35. It is also steeper than the iBaldry fc Glazebrookl (120031 ) value of F = 1.15 ± 0.2 
derived from galaxy luminosity densities in the UV to NIR. It is however well within their 
95% confidence limit of F < 1.7 as well as their measurement of F = 1.2 ± 0.3 based on 
the Ha luminosity density. The difference between their two results suggests that the Ha 
and mid-UV to optical fluxes may have different sensitivities to massive stars. IScald (Il998l ) 
estimated the uncertainty in F, either due to measurement uncertainties, real IMF variations 
or both, in his star cluster IMF based on the spread of results in the literature. Our result 
is well within his range of uncertainty in both mass regimes- F = 1.7 ± 0.5 for 1 — 10 M© 
and F = 1.3 ± 0.5 for 1 - 100 M©. 



6.1. Luminosity Effects 

The luminosity of a galaxy could potentially have an effect on the IMF within it. For 
one the ambient radiation field is likely higher in more luminous galaxies. Figure [6] shows the 
best fitting IMF and x^ values as a function of M^ o.i for all 130,602 galaxies. The galaxies 
have been binned in Mj-o.i such that there are 500 objects in each bin. The bin size was 
chosen to maximize coverage in Mr,o.i yet still keep the random errors in each bin small. The 
solid lines represent the lower and upper 95% confidence region determined from the MC 
simulation. 

Figure [6] reveals a constant value of F ~ 1.4 for galaxies with Mj, o.i between —21 and 
—22 with linear increases in F for both brighter and fainter galaxies. There is also a sudden 
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downturn in T values for galaxies fainter than Mr,o.i = —16.5. Given the sizes of the random 
errors the differences in F between M^^q.! = and —22 are substantial, from 1.59 to 1.41, 
and statistically significant. The agreement with the Salpeter slope is the best for galaxies 
between —21 and —22 in M^^o.!- 

In many ways it is not surprising that previous investigators have not foun d this trend. 
The Milky Way is thought to have a luminosity of My = —20.9 (jDelhaydll965l ): the V and 
tq i filter curves cover roughly the same wavelengths. At comparable luminosities our results 
are similar to Salpeter. The galaxies in the K83 sample have a median Mr = —20 .9 with 
only 16% (18 objects) fainter than M% = -19.7 jEfstathiou. Ellis, fc Peterson! Il988h . This 
is a significant bias toward more luminous galaxies where our results are in ag reement with a 
univer sal IMF. By contrast 30% of our sample is fainter than M*q ^ = —20.44 (IBlanton et al. 
2003bl ). We have a sample of 39,350 galaxies fainter than L*. 



The lower panel of Figure [6] shows that the relative quality of the fits rapidly deteriorates 
as the luminosity of the galaxies decrease. For the brightest galaxies floats around 0.15, 
while in the faintest bin it is over 6. For compariso n My = —18.5 for the Large Mag ellanic 
Cloud and —17.1 for the Small Magellanic Cloud ( ICourteau &: van den Berghlll999l ). This 
trend could indicate that a universal IMF is a good fit to the most luminous galaxies, but 
dwarf galaxies cannot be described by a universal IMF, even if a different universal slope is 
allowed. However it could have a more mundane explanation. It could be that errors are 
over or underestimated as a function of luminosity. It also could arise from deviations from 
our assumption of smoothly varying SFHs. 

We cannot bin our data by stellar mass without assuming an IMF which is contrary to 
the goals of the project. We can repeat the analysis of Figure [6] using M^^o.i in the place of 
Mr,o.i- M^^o.i, being redder, is a better proxy to stellar mass. The resulting plot is nearly 
identical to Figure El which shows that the relationship persists across several wavebands. 

Figure [6] reveals a clear, statistically significant trend in F and with respect to lu- 
minosity. The rest of this section focuses whether this trend is a manifestation of true IMF 
variations or if it is the result of sample biases or poor assumptions. 



6.2. Sources of Bias 

If the IMF is truly universal and our method successfully probes the IMF any subsample 
of galaxies that we could choose should yield the same F value as any other in spite of any 
selection biases or aperture effects. Figure [6] clearly shows that the preceding statement is 
false. In this section we set aside the possibility of IMF variations and search for biases in 
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our sample. 

6.2.1. Magnitude Limited Sample 

Figure El shows that the overall result of F = 1.4437 is really a weighted average. The 
SDSS MGS is a magnitude limited sample, one defined by fiux limits, with both upper and 
lower limits. Table [1] gives the number of objects in each luminosity bin. There are 41,411 
galaxies with -21.5 < Mr,o.i < -20.5 but only 28 for which -14.5 < Mr,o.i < -13.5. As 
such the overall result is heavily biased by more luminous galaxies. 

Malmquist bias will affect any magnitude limited sample. Because brighter objects 
can be seen at greater distances a magnitude limited sample contains bright objects from 
a greater volume of space than fainter objects. The result is that the ratio objects by 
luminosity in a magnitude limited sample differs from the true ratio in nature; brighter 
objects are over-represented. 

To eliminate Malmquist bias volume limited bins where subsamples are complete for a 
range of luminosities are constructed. Figure [7] details the construction of these bins. Given 
both the upper and lower flux limits of the MGS (15.0 < r < 17.77) only a factor of 13 
in luminosity falls in the sample at any given redshift. The redshift limits of each volume 
limited bin are deflned such that no galaxies within the magnitude limits of the bin are 
affected by the flux limits of the MGS. Within each box in Figure [7] the true ratio of galaxy 
luminosities is preserved and is thus free of Malmquist bias. 

Figure [S] shows the results for volume limited magnitude bins. Most error bars are 
smaller than the plotting symbols due to the larger number of galaxies, 329 to 29,701 as 
given in Table [1], in each bin. Figures [6] and [8] show the exact same trends. The largest 
difference in F between the whole and volume limited samples is 0.0116 in the M^^o.! = ~17 
bin. The other notable difference is that the fainter galaxies have larger values in the 
volume limited case. However Malmquist bias across bins is not responsible for the luminosity 
trends in Figure [61 

6.2.2. Redshift 

Another effect of magnitude limited samples is that the faintest galaxies are much closer 
than the most luminous ones. The mean redshifts of the volume limited magnitude bins 
range from z = 0.013 for Mr,o.i = —17 to z = 0.168 for M^^o.i = —23. This corresponds to 
a difference in age of around 1.8 Gyr. As aforementioned, model tracks reveal that age is 
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largely orthogonal to the IMF in our parameter space, but there could be other effects tied 
to age and distance. In addition the IMF could evolve with time. 

The large number of galaxies in our sample affords us the luxury of investigating the 
effects of luminosity and redshift simultaneously to obtain a better understanding of what 
role, if any, the redshift plays in our analysis. Figure [9] shows our fitted parameters for 
all 130,602 galaxies in bins that are 0.25 magnitudes wide in luminosity and 0.005 wide in 
redshift. The upper left panel shows the best fitting F for each two dimensional bin. On the 
upper right the width of the 95% confidence region in F for each bin is shown. At bottom 
left is the log and at bottom right is the log of the number of galaxies in each bin. The 
white contour demarcates the region in which each bin contains at least 50 galaxies. The 
number 50 is arbitrary but it shows the region where Poisson errors are expected to be small. 
The black areas are regions where there are no galaxies with the given parameters. 

Using the plot of F at the upper left we can look for potential redshift biases. This is 
complicated by the fact that at a fixed luminosity there is a limit to the range of redshifts 
in the sample due to the flux limits of the sample described earlier. Looking at vertical 
slices through the plot at any fixed luminosity there is a trend towards larger values of F 
with increasing redshift. However, for horizontal slices of fixed redshift the same relationship 
between F and luminosity that is present for the whole sample is seen modulo a normalization 
factor. 

The right half of Figure M shows a strong relationship between the number of galaxies 
per bin and the width of the 95% confidence region in F. This simply reflects the fact that 
larger samples are less affected by Poisson errors. 

The lower left panel of Figure M provides an excellent example of why our metric of flt 
quality, x^, is so important. Bins with similar numbers of galaxies and ST values can have 
vastly different values of x^- It is worth reminding that the contours in are logarithmic. 
At flxed luminosity the galaxies are better flt by a universal IMF at higher redshift. Similar 
to the sample as a whole the quality of flt improves with luminosity. 

While there does appear to be some weak trending of F and with redshift, redshift 
effects are not driving the relationship seen between IMF and luminosity as it persists at 
flxed redshifts. 
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6.2.3. Aperture Effects 

One explanation for the trend in F with redshift is aperture effects. Again, if the IMF 
is truly universal aperture effects should not exist. The SDSS spectra have a fixed aperture 
of 3" for all galaxies. Depending on the angular extent and distance to a galaxy a different 
fraction of the total light of the galaxy will fall into the aperture. The problem is mitigated 
by the fact that the most distant galaxies are the most luminous and more likely to have a 
larger physical size. As the physical area contained in the aperture increases with distance, 
so too does the size of the galaxies being observed. However, the two effects do not exactly 
balance out. Table [T] shows that the mean aperture fraction for the M^o.i = ~17 bin is 0.20 
and increases to 0.27 at Mr^o.i = —23. On average 35% more of the most luminous galaxies 
fall within the aperture compared to the faintest. 

Figure [10] shows the behavior of our fitted parameters for two dimensional bins of lumi- 
nosity and aperture fraction in the same manner as Figure [9] did for luminosity and redshift. 
For fixed luminosities increasing aperture fraction leads to decreasing values of F. However 
at fixed aperture fraction the qualitative IMF-luminosity trend remains. The values are 
a strong function of luminosity, but does increase with aperture fraction at each fixed 
luminosity. 

The trend with aperture fraction is the exact opposite of what would be expected in 
the presence of a systematic effect operating given the redshift result in Figure [91 The 
nearest galaxies should have the smallest aperture fractions in a particular luminosity bin. 
The nearest galaxies in Figure [9] have the smallest values of F while the smallest aperture 
fractions in Figure [TOl have the largest values of F. 

Figure [10] suggests that the measured IMF is more dependent on the aperture frac- 
tion than the redshift. There are several possible physical explanations for IMF trends 
with the aperture fraction, all of which are related to radial gradients in disk galaxies. 



Padoan. Nordlund. fc Jones! (119971 ) make the theory based claim that the IMF should be a 
function of the original local temperature of the star-forming molecular clou ds. Metallicit y 
gradients are al s o kno wn to exist in disk galaxies, including the Milky Way (jMayor 1976). 



RoUeston et al.l (120001 ) measure a linear, radial light metal (C, O, Mg & Si) abundance gra- 
dient of —0.07 ±0.01 dex kpc~^ in the disk of the Milky Way. Given the increased efficiency 
of cooling with metal lines we would expect the most low mass stars where metallicity is 
the highest- on average towards the center of galaxies. The trend in F in Figure [TD] is 
qualitatively consistent with this idea. 



If there are radial IMF gradients in galaxies one would expect the fits to decrease in 
quality with increasing aperture fraction. A blend of IMFs will not be fit as well as a universal 
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one given our technique. This idea is consistent with the results in Figure [TOl However, in 
well-resolved stellar populations there is no evidence for a relationship betwee n the IMF an d 
metallicity, except perhaps at masses lower than those probed by our method (lKroupall2002l ). 
If metallicity plays a role in determining the IMF the effects are only being revealed as a 
global trend in our large sample of integrated stellar populations. For individual clusters 
metalhcity must play a secondary role to stochastic effects. 



6.2.4- Extinction Correction 

The extinction correction is another potential source of bias. It is possible that there 
is a second order correction that our fairly simple extinction correction fails to take into 
account. This could potentially lead to an erroneous IMF trend with extinction correction. 
This affects the luminosity results because more luminous galaxies tend to be dustier, as 
evidenced in Table [H The problem is further complicated by the fact that dust is thought 
to play an integral part in star formation so it is not unreasonable that an observed IMF 
trend with extinction may be real. 

Figure [11], similar to Figures [9] and [T0|, shows the results of our analysis for two dimen- 
sional bins of luminosity and the extinction correction that was measured and applied. Verti- 
cal slices through the upper left panel of Figure [TT] show that F does depend on Ay,;, trending 
towards lower F values with increased extinction over the region where the Poisson error in 
F is reasonable. Yet again, horizontal cuts of fixed extinction show the IMF-luminosity 
relationship. 

The decreasing F values with increasing extinction is counter-intuitive. Dustier regions 
tend to be more metal rich. If metal cooling plays a significant role in the IMF the dustiest 
regions should have the steepest IMFs. 

At fixed luminosity increases with extinction. As aforementioned our calculated er- 
rors in color and EW (equations [5] and [7]) have a functional dependence on the observed 
emission line extinction. In both cases it is the term proportional to ci/Av,/ which is on 
average the major contributor to the calculated error. Because luminous galaxies tend to be 
more heavily extincted they will also be more likely to have larger errors. This is potentially 
problematic for our observed IMF trend with luminosity. If we are unknowingly underesti- 
mating the errors for faint galaxies with low extinction the source of the poor fit qualities 
of these galaxies could be systematic instead of astrophysical. However the lower left panel 
of Figure [11] shows that the most extincted galaxies have the poorest fits where such a bias 
would suggest that they should fit the best due to the large accommodating errors. 
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6.2.5. Multiple Parameter Biases 

It is also possible that biases in our F measurements could depend on two parameters 
simultaneously. Figure [12] shows the measured galaxy of F as a function of both aperture 
fraction and measured emission line extinction for six volume limited luminosity bins. 

When holding all other parameters fixed, increasing the aperture fraction leads to lower 
values of F in all statistically significant areas of Figure [121 This is the same relationship 
found in the earlier section on aperture fraction. Decreasing values of F are also seen for 
increasing extinction when all other parameters are constant. A notable exception to this is 
that galaxies with large extinction and small aperture fractions favor higher F values. 

Most importantly when looking at a particular combination of aperture fraction and 
extinction the IMF becomes shallower with increasing luminosity until the highest luminosi- 
ties where it becomes steeper again. Even in the narrowest slices of the data set the same 
IMF- luminosity trend is seen, albeit with slightly different absolute values of F. 

6.2.6. Star Formation Strength 

As discussed previously we have allowed two classes of star forming galaxies into our 
sample. 111,806 galaxies (86%) fall in the star forming class and the other 18,796 (14%) 
belong to the low S/N star forming class where the O III or N II lines are weak, but the Ha 
and H/9 lines still have S/N > 5. By comparing the results from these two subsamples we 
can investigate a possible bias of the results with respect to the level of star formation. 

Figure [13] shows the results for both classes as a function of luminosity. As it comprises 
86% of the total sample it is not surprising that the results for the star forming class are 
similar to those of the sample as a whole in Figure [S] 

The low S/N class exhibits a similar qualitative behavior to the set as a whole with 
a few notable differences. The F values are offset by at least 0.08 towards larger F. The 
measured IMF turns toward steeper values at lower luminosities than for the sample as a 
whole. The are several times lower as well. 

While the galaxies in low S/N class meet the same requirement of S/N > 5 in the Balmer 
lines as the star forming class they are biased towards noisier Ha line measurements. This 
corresponds to lines that are either weak (low SFR) or weak compared to the continuum 
(low present SFR compared to the past) both of which lead to low EWs. Another issue 
at play is that the relationship between the Ha line flux and the SFR is dependent on the 
IMF. At a fixed metallicity and SFR increasing F by 0.05 reduces the Ha flux by 20%. In 
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fact the Ha flux of a galaxy with T = 1.00 will be 33 times larger than a galaxy with the 
same SFR and F = 2.00. In the presence of real IMF variations at any fixed luminosity the 
low S/N class will be biased towards galaxies with steeper IMFs. Both low SFRs and steep 
IMFs potentially lead to low S/N Ha flux. However it is difficult to determine the level of 
influence of each effect. 

As shown in figure [1] low Ha EWs lead to larger values of F for any fixed color. As 
the low S/N class tends toward noisier Ha fluxes and therefore EWs it is easy to see from 
equations El and [7] that the errors for this class will tend to be larger. This in turn leads 
to lower values. 

As the qualitative IMF-luminosity trend occurs in both star forming classes the strength 
of star formation is unlikely to be a significant bias on our results. 



6.2.7. The f Ratio 

As mentioned in the data corrections section, the / ratio is ratio of the extinction 
experienced by the nebular emission lines to that experienced by the stellar continuum. The 
assumption of a value for / could potentially bias our results. An alternative way of looking 
at the same problem is that our data in the color-Ha EW plane can be used to constrain 
the / ratio by assuming a universal Salpeter IMF. 

Figure [H] gives the results of this analysis for the data set as a whole. A value of / = 2.0 



i s foun d with x'^ = 1-068. This is in good agreement with to the lCalzetti. Kinney. &: Storchi-Bergmann 



(119941 ) value of / = 2.0 + 0.6/ — 0.4. The quality of the fit in the best case is worse than 
in Figure [51 Part of the reason for this is that the errors used were slightly smaller as the 
0"/ terms in equations [5] and [7] are set equal to 0. The quality of the fit drops sharply below 
/ = 2 and more gradually for / > 2. Values of / near 1 are heavily rejected. However in 
this particular plot the results are dominated by luminous galaxies. 

The values of / as a function of luminosity are shown in Figure [T5l For galaxies M^ o.i = 
— 19 and brighter the best value of / is consistent with / = 2. The faintest two bins the 
prefer an / ratio closer to 2.5. However the lower panel shows that this new / value does 
not translate to improved fit quality. In fact the faintest galaxies have in general smaller 
measured extinctions and are therefore less susceptible to changes in /. The same qualitative 
trend of worsening fits with decreasing luminosity seen when allowing F to float is seen with 
a varying / value. 



Together these two / ratio plots provide a number of insights. For one it shows that our 
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choice of the / ratio is very sensible and provides an independent confirmation of other / ratio 
measurements. The fact that our best fitting F values are at least 0.05 above the Salpeter 
value cannot be reconciled by changing the geometry of the dust screen. It provides further 
evidence that the relationship between F and luminosity is not a function of extinction or a 
byproduct of our extinction correction. 

6.2.8. Summary 

In this section we have investigated several possible sources of bias to account for our 
observed trend between the IMF and luminosity. Relationships between the IMF and red- 
shift, aperture fraction, extinction and star formation strength have been uncovered. Two 
parameter biases were also found. 

In all cases in narrow slices through the data where potential biases are held fixed the 
qualitative IMF-luminosity relationship appears. The parameters primarily act to offset the 
value of F at a particular luminosity. The ratio of continuum to emission line extinction, /, 
was found to be a sensible choice and the results are not sensitive to small changes in this 
value. 

There are two possible interpretations to the relationships between the IMF and poten- 
tial biases. One is that they are systematic effects due to some problem with our measurement 
of F. The second is that they are real physical effects. It is not clear from the data which of 
these statements is more correct. 



6.3. Star Formation History 

In the previous section several possible sources of bias were investigated, but none were 
able to account for our observed trend in F with luminosity or the inability of a universal 
IMF to fit low luminosity galaxies. Figure [16] shows the distribution in color-Ha EW space 
for the least and most luminous bins, M^ o.i = and —23. From this figure it is apparent 
that the most luminous galaxies lie roughly parallel to the IMF tracks while the faintest 
galaxies are more perpendicular to the tracks. In the low luminosity bin there are galaxies 
which are simultaneously blue and have low EWs. These galaxies are not consistent with a 
universal IMF with F = 1.35 and as mentioned before are not consistent with a universal 
IMF with a different slope. In addition the faintest galaxies have the lowest extinctions they 
are the least sensitive to dust and / ratio issues. Before concluding that this is evidence 
for IMF variations we must first consider whether our model assumption of smooth SFHs is 
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justified. 

6.3.1. Effects of Star Formation Bursts 

The SFH of individual galaxies is tlie most problematic aspect of tfie K83 analysis. A 
sudden burst on top of a smoothly varying background will immediately increase the Ha 
EW. This is due to the formation of O and B stars which indirectly increase the Ha flux 
through processing of their ionizing photons. The new presence of O and B stars also makes 
the color of the galaxy bluer. Both of these effects are proportional to the size of the burst. 

After the burst is over the Ha EW is smaller and the colors are redder than they would 
be if the burst had not occurred. The Ha EW drops because there is no longer an excess 
of O and B stars and their ionizing photons which reduces the Ha flux to pre-burst levels. 
However there is now an excess of red giants due to the less massive stars from the burst 
leaving the main sequence. This increases the continuum around the Ha line which further 
drops the EW in addition to making the galaxy colors redder. After enough time has elapsed 
after the burst the galaxy returns to the same position in the color-Ha EW plane it would 
have occupied had no burst occurred, although it will have taken longer to get there. 

Figure [T71 gives one example of this cycle. A solar metallicity galaxy with F = 1.35 and 
an exponentially decreasing SFH with r = 2.15 Gyr experiences a burst of star formation at 
an age of 4.113 Gyr which lasts 250 Myr and forms 10% of the stellar mass. The black dots, 
spaced at 100 Myr intervals, show that comparatively more time is spent below the nominal 
track than above it. The peak Ha EW is reached just 5 Myr after the start of the burst. If 
you happen to be observing the galaxy during the burst a shallower IMF will be measured 
(assuming a burst-free SFH), after the burst for 1 Gyr a steeper IMF will be measured and 
after that the effects of the burst largely disappear, although the galaxy will appear younger 
than it actually is. 

Figure [TS] gives the best fitting IMF slope from our analysis, which assumes no bursts, 
as a function of age for the galaxy in Figure [T71 The jitter in the best fit F values is due to 
the discrete nature of our model grid, and the fact that the model tracks for different IMFs 
run together at large ages. Within a 300 Myr period during and just after the burst the best 
fitting IMF slope is anywhere from F = 1.00 to 1.95. This shows that even if a universal 
IMF exists the SFH can mimic a huge range of IMF models. Roughly 1 Gyr after the burst 
the measured IMF is back to its true value. A galaxy with a bursty SFH viewed at a random 
time will be biased toward a steeper IMF than what is the true IMF. While one model is 
not an exhaustive study of the effects of the SFH on IMF measurements it does give a good 
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sense of what issues arise. 

To eliminate this uncertainty we investigated cutting the sample on SFH. In order to 
detect a relative lack or excess of present star formation it is necessary to measure both 
the present and past star formation rates or at least be able to compare the two in some 
way. The problem lies in the fact that conversions of observables into star formation rates 
assume an IMF to do so. Our aim is to measure the IMF so we cannot make strong a priori 
assumptions about it. Instead of biasing or results from the start we fit all galaxies and then 
try to determine the affect of SFHs on our conclusions. 



6.3.2. Single Burst Models 

The simplest burst model is that of a single burst at a random time on top of our 
smooth exponential SFHs. A grid of 1000 SFHs was constructed by first selecting one of the 
24 smoothly varying SFHs at random. A burst lasting 200 Myr was superposed on the SFH 
at a time selected uniformly at random over a range of 12.5 Gyr. The strength of the burst 
was randomly selected up to 40% of the total stellar mass, with preference given to smaller 
bursts. The colors and EWs of these SFHs were calculated at 1 Myr intervals over 12.5 Gyr 
for the IMFs F = 1.35 and 1.80 and a fixed metallic ity of Z = 0.01. This metallicity choice 



is based on the luminosity-metallicity relationship in lTremonti et al.l (l2004j ) for Mr o.i = —17 
galaxies. Plotting these models in the color-Ho; EW plane shows that all observed data 
points are covered by either IMF. 

To test whether the observed distribution of points in the M^o.i = — 17 bin can be 
explained by bursting SFHs the MC techniques of §5 were used. The models of §4 were 
replaced with the grid of single burst models. The 100 MC simulations were constructed 
as described earlier, but using only the 329 Mr,o.i = —17 galaxies as a basis. The analysis 
yielded = 0.027 for F = 1.35 and = 0.020 for F = 1.80, both of which are over fits. 

This shows than an individual galaxy can be fit with an arbitrary IMF given the freedom 
to choose a SFH. However, our advantage is that we have many galaxies and the distribution 
of the properties of the best fit models can be shown to be implausible. 

Figure [I9] demonstrates the problem with the single burst model. On the left of the 
figure the distribution of the best fitting burst strengths are plotted as a fraction of the total 
stellar mass formed. At right is the distribution of the best fitting times from the burst 
onset. For example, a galaxy which best fits a model with a burst at 1.000 Gyr at an age 
of 1.211 Gyr has a time from burst onset of 211 Myr. This measure is used because in an 
investigation of the effects of bursts the age relative to the burst is more important than the 
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age given that the bursts occur at different, random times across the models. 

For both F = f .35 and F = 1.80 the number of objects best fit by a model prior to the 
burst is 11 + 3/ — 5, or 3%, and are not plotted. In both cases the distributions of the best 
fitting ages and ages at which the burst begins are roughly uniform. Given this fact it is 
expected that half of the galaxies should be best fit by a pre-burst model. Furthermore the 
right hand panel shows sharp, significant discontinuities in the distribution of best fitting 
time from burst onset. Again, viewed at random times this distribution should be uniform 
but is highly peaked in the 25 Myr at the start of the burst and the 25 Myr just after the 
burst ends. In both cases the errors bars show that the discontinuities are significant. In the 
case of F = 1.35, 3.5 times as many galaxies are in the 200 Myr after the burst ends than the 
200 Myr during it and this 400 Myr accounts for 57% of all galaxies. Although our sample 
is r-band selected the stars in the 0.7 to 3 Mq range which dominate the red continuum in 
the red giant phase do not start to leave the main sequence for 300 Myr. The sharp increase 
in galaxies fit at 200 Myr after the burst cannot be due to a selection effect. Assuming a 
universal IMF this points to a strong coordination of SFHs across a population of galaxies 
unrelated in space. These arguments show that while a single burst model can fit the data 
extremely well, it does not do so in a physically self-consistent fashion. 



6.3.3. Multiple Burst Models 



To find a physical motivation for SFH models for low luminosity galaxies we look to- 
wards the Local Group. There have been a number of recent studies of the SFHs of local 
dwarf galaxies which use HST to get color- magnitude diagrams (CMD) of resolved stellar 
populations. The SFH is determined by fitting isochrones to the CMD. The sample here is 
biased by Local Group membership and by what galaxies have been observed to date. The 
galaxies mentioned here give a point of reference rather than a well-defined distribution of 
SFHs. 

The blue compact dwarf (BCD) UGCA 290 was found to quiescently form stars over 
the past Gyr up until a ten-fold increase in SF R from 15 to 10 My r ago which more recently 
has decreased to a quarter of its peak value (ICrone et al.ll2002l ). The dwarf irregular IC 
1613, which is relatively isolated and non-interacting, was found to have SFR enhance d by a 
factor of 3 from 3 to 6 Gyr ago without evidence of strong bursts (jSkillman et al.ll2003l ). The 
dwarf irregula r NGC 6822, also relatively isolated, is found to have a roughly constant SFH 
(IWyderl |2003| ) . The BCD NGC 1705 is found to be gasping- a SFH marked by mo derate 
activity punctuated by short periods of decreased star formation (lAnnibali et al.ll2003l ). The 
authors also note that NGC 1705 is best fit by an IMF with F = 1.6. NGC 1569 likely 
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experienc ed three strong bursts in the last Gyr as well as a quiescent phase from 150 to 300 



Myr ago (lAngeretti et al.ll2005l ) 



Inspired by the preceding Local Group SFHs we constructed six SFH classes with multi- 
ple bursts. These SFH classes are described in Table El Each class starts with an underlying 
smooth SFH. SFHs #1 & 2 have no star formation, #3 & 4 have a constant SFH and ^^5 
& 6 have exponentially decreasing SFHs like those previously described. Star formation 
discontinuities are then superimposed on top of the smooth SFHs. These discontinuities are 
in the form of increased (bursts) or decreased (gasps) star formation for periods of 200 Myr. 
The time and spacing of the discontinuities is random with the mean interval between bursts 
listed in Table [2l 

For each SFH class described above we randomly generated 1,000 SFHs. Colors and 
EW widths were calculated for each SFH using F = 1.35 and Z = 0.01. According to the 



SDSS mass-metallicity relationship (ITremonti et al.ll2004l ) galaxies at Mr,o.i = —17 will on 
average have Z = 0.01. 

We then repeated our analysis with the 100 MC simulations of the Mj.,o.i = — 17 
galaxies in the same manner as for the single burst models. The results of our analysis for 
each of the six SFH classes are shown in Figure [201 SFH 5, the gasps on top of exponential 
SFHs, is the best fit with = 0.09. Extended periods of no star formation punctuated by 
bursts (SFHs # 1 & 2) do not fit the data. 

As was the case for the single burst models an unreasonable fraction of galaxies are 
best fit by SFHs in the 20 Myr immediately following a burst or the first 20 Myr of a gasp. 
If the SFH models are reasonable we should see roughly equal numbers of galaxies in each 
time bin. There is no reason why all of the low luminosity galaxies across the large volume 
of space in the SDSS footprint should have experienced coordinated bursts. However each 
panel of Figure [201 has at least 40% of the galaxies in one 20 Myr bin. 

One explanation for this is that it is an artifact of our sample being selected in the r 
band. However spectral synthesis models show that for instantaneous bursts of star formation 
the ro.i magnitude is brightest at the burst time and decays smoothly for a range of F. If 
anything it is more likely to catch galaxies during a burst rather than after or after a gasp 
instead of during one. 

Regardless of the SFH model the presence of blue galaxies with low Ha EW requires 
a recent discontinuity in the SFR for F = 1.35. Based on the evolution of the ro.i band 
luminosity we expect to see a similar number of galaxies with excess Ha EWs. The fact 
that these galaxies are missing shows that the discontinuous SFH models do not match 
our observations. Therefore IMF variations are a more likely explanation for the observed 



-35- 



distribution of Mr,o.i — —17 galaxies. 

6.3.4- Recovering T from Synthetic Data 

As a last exercise the best fitting SFH models from the previous section can be run 
forward to see if the correct IMF can be recovered. For each SFH model grids 10,000 data 
points were chosen by selecting a random SFH and a uniformly distributed age. Normal 
errors were added using the error characteristics of the Mr,o.i = —17 bin. This synthetic 

data was analyzed in the same fashion as the real data in the earlier sections. For the 
single burst models the recovered IMF models for F = 1.35 and 1.80 were 1.34 and 1.79 
respectively with best fitting = 0.80 and 0.65. For SFH 5 the recovered IMF was also 
1.34 with = 1-5. In all three cases the correct IMFs were recovered although the fit was 
worsened by by the burst activity. This reinforces the difficulty in producing enough blue, 
low Ha EW galaxies to match the observed data with simple SFH models. 

7. H5a Absorption 

In the previous sections we have expanded on the K83 method and exploited the Ha 
and color information as much as possible. In the bias section we found that various possible 
biases do not fully explain either the increased values of F or the poor fit to a single IMF in 
the lowest luminosity bins. In the SFH section we found that an arbitrary F value coupled 
with a plausible SFH with bursts or gasps can account for the position of any individual 
galaxy in the color-Ha EW plane. However taking the population of Mr,o.i = — 17 galaxies 
together necessitates an incredibly unlikely coordination of SFHs across the disparate group 
of objects. This points to the extraordinary conclusion that while the IMF may be universal 
across luminous galaxies, it is not in fact universal in low luminosity galaxies. Such an 
extraordinary claim would ideally be backed by extraordinary evidence. In this section we 
take a look beyond the K83 method for some reinforcement of our result. 

The B.6 absorption feature can be used to gain additional insight into the nature of 
stellar populations. H5 absorption is due to absorption lines form stellar photospheres. The 
Balmer absorption lines are most prominent in A stars and weaken due to the Saha equation 
for both hotter and cooler stars. As such the B.6 absorption is a proxy for the fraction of light 
of a stellar population being supplied by A stars, and to a lesser extent B and F type stars. 
In a stellar population of a uniform age the R6 absorption will peak after the O and B stars 
burn out, but before the A stars leave the main sequence. For this reason the strength of 
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US can be used help determine the age of a population or to detect bursts of star formation 
which occurred around 1 Gyr in the past. 



Worthey &: Ottavianil (119971 ) describes two different methods for measuring H5 absorp- 
tion. The H^F definition is tuned to most accurately measure the H5 absorption from F 
stars. The H(5a definition has a wider central bandpass to match the line profiles of A stars. 
They state that the H5a definition is better to use for galaxies because it is less noisy in low 
S/N galaxies and velocity dispersion acts to widen absorption features. On the downside 
the narrower H(5f definition is much more sensitive to population age where it can be used. 



Worthey fc Ottavianil (119971 ) observationally determines the range of H5a values to be from 



13 for A4 dwarf stars to -9 for M-type giant stars. 

H5a values c an be measured from the SDSS spectra. Like the Ha EW values our H(5a 
values come from Tremonti et al. J2OO4I ) instead of the SDSS pipeh 



nne. 



Figure [21] compares the distribution of the B.6a values for the M^^o.! = —17 and Mr,o.i = 
—23 luminosity bins. For reference recall that Figure [16] plots the color and Ha EWs for 
these two bins. The difference between the high and low luminosity bins is clear. Gaussian 
profiles can be fit to both distributions. For the Mr,o.i = —23 bin the standard deviation 
of the profile is twice as large as the measurement error in H5a suggesting that the true 
distribution has a range of values. Assuming both the errors and underlying distribution are 
Gaussian the distribution of H5a for the luminous galaxies is centered at H(5a = 3.4 with 
a standard deviation of 1.5. By contrast, the Gaussian fit to the M^^o-i = —17 bin has the 
same standard deviation as the median error in the galaxies. This is consistent with nearly 
all of the galaxies having B.6a = 6.1. 

This shows that the fainter galaxies on average have significantly larger fractions of A- 
type stars amongst their stellar populations. It also shows a seemingly unlikely coordination 
of H5a in the low luminosity galaxies, reminiscent of the earlier suggestion of coordinated 
SFHs. The question then becomes why? 

For a possible explanation we look again to the models. Modeling the behavior of H(5a 
requires an extra step. The standard PEGASE.2 models do not have the required resolution 
to ac curately measure the H ^a index. This is remedied with the use of the PEGASE-HR 



code (iLe Borgne et al.|[2004l ) which uses a hbrary of echelle spectra of 1503 stars to calculate 
spectral synthesis models with R=10,000 over the range of 4000 to 68OOA. Using PEGASE- 
HR in the low resolution mode yields the same results as PEGASE.2, and the same input 
parameters are used for both codes. The models here are the same as those described 
earlier in the models section, but have been recalculated using PEGASE-HR to allow B.6a 
measurements. 
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Figure [22] shows the behavior of the H5a index as a function of age for four different 
IMF models, F = 1.00, 1.40, 1.70 and 2.00. For each IMF models of all metaUicities and 
smooth SFHs are plotted for each age. The qualitative behavior of all models is the same. 
The H5a holds steady for the first 20 to 40 Myr before increasing to a peak value at 700 Myr 
to 1 Gyr and then falls off. Prior to reaching the peak value for each IMF the metallicity 
has the strongest effect on the B.6a index. At this point the lowest metallicity galaxies have 
the highest US a- After the peak the SFH has the strongest effect with the constant and 
increasing SFHs maintaining higher USa values. 

After a few Gyr the differences in B.6a between models with different IMFs disappear. 
Prior to that there are three main differences. First, the peak B.6a values are higher for larger 
values of F. For F = 2.00 H5a reaches a maximum value of nearly 8 and the maximum value is 
similar for all metaUicities. For F = 1.00 the maximum value ranges from 4.5 to 6 depending 
on the metallicity. This is due to the fact that the steeper IMFs have fewer luminous massive 
stars to dilute the B.6 absorption features from the A star population. Secondly, the low F 
models have H^a values that start their initial increases at a later time. Lastly the low F 
models reach their peak values later in time than those with fewer massive stars. 

To compare the B.6 a values for the Mr,o.i = — 23 and —17 bins to the models the range of 
the middle 90% values from Figure [21] for each bin are overlaid on Figure [221 Once again the 
Mr,o.i = —23 bin is in good agreement with our assumption of a universal IMF and smooth 
SFH. The range of H5a values can be accomplished with a single Salpeter-like F = 1.40 IMF 
with only the proviso that most galaxies be older than a few Gyr or younger than 300 Myr. 
However the exact same statement can be made for 1.0 < F < 2.0 so H(5a provides only a 
constraint on the age of the most luminous galaxies, but not the IMF. 

For the M^^o-i = ~17 galaxies the distribution of B.6a cannot be achieved with the 
shallowest IMFs investigated under the assumption of smooth SFHs. However Salpeter and 
steeper IMFs can be accomplished. What changes is the range of ages over which the models 
have the correct B.6 a- Steeper IMFs require that the galaxies be either older than a few Gyr 
or 200 Myr old to accommodate the observed B6a values. Most troubling is that there does 
not appear to be any reason why B6a should stack up at 6.1 for the M^o.i = ~17 galaxies. 

H^A values can also be calculated for the same SFH class models used earlier. Unfortu- 
nately this does not provide any added constraints. B6a values remain elevated for several 
100 Myr after the start of a burst or the end of a gasp of SFH. The behavior of B6a in the 
presence of SFH discontinuities provides no need for galaxies to stack up in the narrow 20 
Myr intervals seen in the earlier multiple burst section, nor do these models suggest why the 
low luminosity galaxies are consistent with a single value of B6a- 
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There is no satisfactory model to account for the H5a distribution of the low luminosity 
galaxies. However the SFH results from the previous section strongly suggest that the 
incredible coordination of discontinuities is highly unlikely. Steeper IMFs do allow for higher 
H5a values for longer periods of time, thus relaxing the SFH coordination requirement. This 
agrees with our earlier results for faint galaxies. The low luminosity galaxies are most likely 
the result of a mix of IMFs which are on average steeper than Salpeter. 

8. Conclusions 

The goal of this paper was to revisit the K83 method for inferring the IMF from inte- 
grated stellar populations and to harness the richness of the SDSS data, improved spectral 
synthesis models and greater computational power available today to make a state- of- the- art 
measurement of the IMF. The quality of the SDSS spectroscopy allowed us to address several 
of the limitations of K83 and KTC94- we resolve the [N II] lines (a significant improvement 
in the accuracy of the Ha EWs of individual galaxies), ehminate contamination from AGN, 
make extinction corrections for individual galaxies and fit underlying stellar absorption of 
Ha. We succeeded in achieving more accurate EWs for individual galaxies. The median 
total EW error for our sample is 17% compared to a 10% uncertainty in EWs combined with 
a 20-30% uncertainty in the extinction correction for K83. 

We expanded the grid of models to allow for a range of ages and metaUicities. We used 
minimization to go beyond differentiating between two or three IMF models to actual 
fitting for the best IMF slope. The vast size of the SDSS sample allowed us to both drive 
down random errors and to cut the data into narrow parameter ranges which were still 
statistically viable. 

The size of the DR4 sample yielded AF = 0.0011 95% confidence region due to random 
error for the sample as a whole. Even the volume limited M,. o.i = —17 luminosity bin with 
only 329 objects has a random error of AF = 0.0086. Only in bins with fewer than 10 
objects do the random errors become significant. Our IMF fitting is therefore dominated by 
systematics. 

Originally we believed that our systematics would be dominated by the efi^ccts of SFH 
discontinuities. However we conducted several experiments where we selected populations 
of galaxies from models with bursting or gasping SFHs and gave them measurement errors 
consistent with those in the Mj.,o.i = —17 luminosity bin. To our surprise our minimization 
revealed the true IMFs with AF = 0.01. The main effect was to reduce the quality of the 
fits. This is due to the fact that Ha EWs return to nominal levels in a relatively short time 
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after SFH discontinuities. 

Another way to estimate the size of the systematic errors is to look at the trends of 
the Mrfi,i = —21 and —22 luminosity bins, because they have the largest membership, in 
Figures [9| [TO] and [TTl Assuming that the IMF is universal and that our method is perfect we 
should get the same answer for any subset of the data we might choose. The largest ranges 
are AF ~ 0.12 for redshift binning, AF ~ 0.19 by aperture and AF ~ 0.19 by extinction. 
Conservatively then the systematic error is ±0.1. 

There are two points to be kept in mind about this estimate of the systematic error. 
For one it is the systematic error in the exact value of F. Even in Figure [12] where more 
narrow bands of measured extinction and aperture fraction are considered the same trends 
with luminosity are seen as with the sample as a whole. The relative systematics between 
luminosity bins in these narrow slices is much smaller. The second thing to remember is 
the way in which we empirically defined our systematic error discounts the possibility of 
IMF variations. What we have called systematics could actually be science. If galaxies have 
radial IMF gradients or if dust content plays a strong role in star formation the systematic 
error could be much smaller. The main area in which we were unable to improve upon the 
K83 and KTC94 studies is that they were able to match the aperture size to the galaxies 
which avoids the issue of aperture effects. 

In spite of a more quantitative approach, like K83 and KTC94 the results are mostly 
qualitative. However there are four key results from our investigation. 

First, for galaxies brighter than M.r,o.i ~ —20 the best fitting IMFs are Salpeter-like 
(F ~ 1.4). In addition the assumption of a universal IMF and smoothly varying SFHs is a 
good fit. This is reassuring as it follows the conventional wisdom and provides confidence 
that the method works. 

Secondly, galaxies fainter than Mr,o.i ~ —20 are best fit by steeper IMFs with larger 
fractions of low mass stars. For these galaxies a universal IMF and smooth SFH is a poor 
assumption. This result is in qualitative agreement with evidence that LSBs have bottom- 



heavy IMFs fiLee et al.ll2004j ) 



Thirdly, while breaking the IMF-SFH degeneracy for individual galaxies using the Ha 
EW and color is hopeless, for a statistical sample of galaxies the degeneracy can be broken. 

Lastly, given our analysis of discontinuous SFHs it appears that the IMF is not universal 
in low luminosity galaxies and fewer massive stars are being created in these galaxies. 

It is worth mentioning the main caveat of our F values again. As illustrated in Figure 
[5] IMF parameterizations are themselves degenerate in our parameter space. Increasing 
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the IMF slope has a similar effect to lowering the highest mass stars that are formed or 
increasing the fraction of intermediate mass stars. This method cannot explicitly determine 
if two pop ulatio n s hav e the same underlying IMF. Figure shows that for the sample as a 
whole the IScald (Il998l ) three part power law yields nearly the same result as our two part 
power law. However our method is sensitive in many cases if the IMFs are different. 

In terms of star forming cloud temperatures the harsher ambient radiation and larger 
number of sources of cosmic rays present in more luminous galaxies agree qualitatively with 
our results. With the extra energy hitting the star forming clouds larger masses may be 
needed for co ntraction and f ractionization may end so o ner, s uppressing the formation of less 
massive stars (lLarsonlll998l ). ICedres. Cepa. &: Tomital (120051 ) find that while the H II regions 
of the luminous grand design spiral NGC 5457 (M 31) can be reproduced by a single Salpeter 
IMF, for the low luminosity flocculent galaxy NGC 4395 a blend of two IMFs is r equired. 
Howe ver, such trends are not seen in studies of well-resolved stellar populations (IKroupa 
2002h . 



Another explanation for the absence of massive stars is that the massive stars are there, 
but are not visible. Extinction to the center of star forming regions, where massive stars 
preferentially exist, can reach Ay ~ 20 (lEngelbracht et al.lll998l ). However the low luminosity 
galaxies have the lowest observed extinctions (see Table [1]) which is the opposite of what 
would be expected given our IMF results. 

It is also possible that the IMF is in fact universal, but the way in which it is sampled in 
embedded star clusters lea ds to an integrated galaxial IMF which varies from the true IMF. 
Weidner fc Kroupal (120051 ) use a universal IMF with the assumption that stars are born in 
clusters where the maximum cluster mass is related to the star formation rate. For a range of 
models this leads to a narrow range for the apparent IMF in high mass galaxies. For low mass 
galaxies the IMF is steeper with a wider range of slopes. The re sults here are in qua li tative 
agreement for some of the integrated galactic IMF scenarios in IWeidner fc Kroupal (120051 ) 
given that there should be a rough correlation between galaxy luminosity and mass. Once 
again, lElmegreenl (120061 ) argues that the galaxy wide IMF should be the same as the IMF 



in clusters regardless. 



In light of the theory of IWeidner fc Kroupal (120051 ). whether the results of this paper 
speak to a relationship between environment and the formation of individual stars is open to 
interpretation. However the impact on the modeling and interpretation of the properties of 
galaxies is clear. iKoppen. Weidner. fc Kroupal (120071 ) note that the integrated galaxial IMF 
is the correct IMF to use when studying global properties of galaxies. Even if the IMF of stars 
is in truth universal it may currently be misused in the modeling of galaxies. Furthermore 
a varying integrated galaxial IMF could open the door to new insights in galaxy evolution. 
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For instance, iKoppen. Weidner. &: Kroupal (120071 ) suggest that the observed mass-metaUicity 
relationship in galaxies naturally arises from a variable integrated galaxial IMF similar to 
the results of this paper. 

Future work will expand in several directions. The success constraining the IMF with 
only two observed parameters (albeit carefully chosen to be orthogonal to systematic errors) 
motivates a more expansive analysis with more parameters. Information from wavelength 
regimes beyond the SDSS can be used. For instance the absorption strength and P-Cygni 
profile shape of th e C IV A1550 line due to massive sta rs is sensitive to the IMF slope and 
upper mass cutoff (ILeitherer. Robert. &: Heckmanlll995l ). However it could be contaminated 



by ab sorption from the interstellar medium and would need to be disentangled (IShapley et al. 



20031 ). A full Markov Chain Monte Carlo analysis could be implemented fitting to multiple 
spectral features and marginalizing over a range of SFHs. In addition to luminosity, surface 
brightness and gas phase metallicity can be tested for systematic IMF variations. While the 
results for luminous galaxies are already dominated by systematics, the continued progress of 
the SDSS can provide better statistics for a more detailed analysis of what physical processes 
are behind the IMF variations in faint galaxies. 
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Table 1. Luminosity Bin Details 
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Note. — Properties of the M,._o.i luminosity bins. Columns 2 through 9 give values for volume limited 
magnitude bins while columns 10, 11 and 12 give values for all sample galaxies within the luminosity range, 
ap is the mean aperture fraction. * The sample as a whole is not volume limited but the more detailed 
information is given for reference. 



Table 2. Multi-burst SFH Models 



Name 


Type 


Length 


Relative 


Spacing 


Underlying 






(Myr) 


Strength 


(Gyr) 


SFH 


SFH 1 


burst 


200 




3 


none 


SFH 2 


burst 


200 




1 


none 


SFH 3 


burst 


200 


4.0 


3 


constant 


SFH 4 


gasp 


200 


0.0 


1.5 


constant 


SFH 5 


gasp 


200 


0.1 


1.5 


exponential 


SFH 6 


burst 


200 


5.0 


1.5 


exponential 



0.0 0.5 1.0 



Fig. 1. — Distribution of all 130,602 program galaxies in the ((7 — ^)0.1 vs. log(Ha EW) plane. 
The contours are logarithmic. Outside the last contour individual points are plotted. The 
thick lines are model tracks with exponentially decreasing SFHs with r = 1.1 Gyr and solar 
metallicity. The age increases along the tracks from 100 Myr in the upper left to 13 Gyr in 
the lower right. The upper line has F = 1.00, the middle line is similar to Salpeter's IMF 
with F = 1.35 and the lower line has F = 2.00. They are identical to the solid lines in Figure 
131 The cross in the lower left indicates the median error bars of the sample. The arrows are 
dust vectors are for typical observed Balmer decrements (Ha/H/? = 4) for different values of 
/. The values of / are 1, 2 and 4 clockwise from the one pointing left. 
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Fig. 2. — Examples of the degeneracy of IMF parameters. All models shown have solar 
metallicity, age 6 Gyr and the tracks are lines of varying SFH from a rapidly decreasing SFR 
on the lower right to slowly increasing on the upper left. In both panels the bold track has 
r = 1.35 and the dashed lines have IMFs F = 1.00, 1.20, 1.40, 1.60, 1.80 and 2.00 from top 
to bottom, all with an upper mass cutoff of 120 Mq. (a) The solid lines are models where the 
upper mass cutoff of the IMF is reduced from 120 Mq to 90, 60 and 30 Mq as labeled, (b) 
The nominal two-piece IMF is replaced by a three piece IMF in the solid hnes. In all solid 
tracks F = 1.35 above 10 M© and F = 0.50 below 0.5 M©, but is altered in the intermediate 
mass region from 0.5 to 10 Mq. The sohd hnes have Fj = 0.50, 0.85 and 1.15 over 0.5 to 10 
Mq as labeled. 
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Fig. 3. — Effects of model parameters in the ((? — t)o.i vs. log(HQ; EW) plane. In both panels 
ages of the stellar population increase along the tracks from 100 Myr on the upper left to 
13 Gyr on the lower right. Also in both panels the upper set of tracks have F = 1.00 while 
the lower set have F = 2.00. The arrows are dust vectors are for typical observed Balmer 
decrements (Ha/H/3 = 4) for different values of /. The values of / are 1, 2 and 4 clockwise 
from the one pointing left, (a) The effect of smooth SFH variation is shown in the left panel. 
The solid lines have exponentially decreasing SFRs with r — 1.1 Gyr, the dotted hnes have 
more slowly falling SFRs with r = 2.75 Gyr and the dashed lines are increasing oc 1 — e~*/^ 
where r = 1.5 Gyr. These tracks have solar metallicity. (b) The effects of metallicity at 
fixed SFH with r = 1.1 Gyr. Z = 0.005 for the dashed line, 0.010 for the dotted hue, 0.020 
for the solid line and 0.025 for the dot-dashed line. The solid lines are identical across the 
panels. This figure demonstrates that the effects of model parameters are largely orthogonal 
to IMF variations. 




- 50 - 




0.0 0.5 1.0 



Fig. 4. — All model points for F = 1.00 (cyan), F = 1.35 (black) and F = 2.00 (magenta). 
Each IMF has 18,480 calculated model values. Z ranges from 0.005 to 0.025, ages range 
from 100 Myr to 13 Gyr and SFHs cover the range described in the text. 
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Fig. 5. — Value of as a function of F for the entire sample (solid line). The "X" marks the 
minimu m value of 6 1671.4 at F = 1.4525. Also plotted are the f or sonae "cla ssic" IMF 
models. ISalpetei Jl955h ffiUed ci rcle) , iMiller fc Scak] Jl979h (asterisk) , IScalol Jl986h (square), 
Kroupa. Tout, fc Gilmore ( 1993 ) (diamond), Scalol ( 19981 ) (triangle) and Kroupa ( 2001 ) (plus 



sign) plotted at rough estimates for equivalent F values using our parameterization. 
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Fig. 6. — MC simulation results for the full sample binned by M^^o-i- Each diamond represents 
500 galaxies plotted at the mean M^^o.i value of the bin. The points are the most dense around 
Mr,o.i = —21.5 as the sample is dominated by galaxies in this luminosity regime. The solid 
lines represent the upper and lower 95% confidence region measured for each bin. Top panel: 
The best fitting F values as a function of r-band luminosity. Lower panel: The values for 
each luminosity plotted on a log scale. 
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Fig. 7. — Definition of volume limited magnitude bins. The dotted lines show the luminosities 
corresponding to the flux limits of the SDSS MGS (15.0 < r < 17.77) as a function of redshift. 
The solid boxes are the volume limited magnitude bins used in this paper. Within each box 
no galaxies within the magnitude range of the box are affected by the flux limits of the 
sample. 
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Fig. 8. — MC simulation results as a function of Mr,o.i where the magnitude bins are volume 
limited. The error bars represent the 95% confidence interval from the simulation and are 
in most cases smaller than the plotting symbol. Top panel: The best fitting F values as a 
function of r-band luminosity. Lower panel: The values for each luminosity plotted on a 
log scale. 
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Fig. 9. — Values of fitted parameters as a function of both luminosity and redsliift. The bins 
extend 0.25 magnitudes in luminosity and 0.005 in redshift. Clockwise from upper left the 
frames show the best fitting F, the width of the 95% confidence region in F from the MC 
simulation, the log of the number of galaxies in each two dimensional bin and the log of y^. 
The shading levels for each panel are given by the adjacent vertical color bars. Black areas 
indicate regions where there are no galaxies with the respective combination of redshift and 
luminosity. The red contour indicates the region in which there are at least 50 galaxies in 
each 2-D bin. 
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Fig. 10. — Values of fitted parameters for 2-D bins cut on luminosity and aperture fraction. 
The bins extend 0.25 magnitudes in luminosity and 1% in aperture fraction. The description 
is identical to Figure [9] as are the shading levels. 




Fig. 11. — Values of fitted parameters for 2-D bins cut on luminosity and emission line 
extinction. The bins extend 0.25 magnitudes in luminosity and 0.06 magnitudes in extinction. 
The description is identical to Figures [9] and [TOl the shading levels. 
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Fig. 12. — Best fitting T values as a function of aperture fraction and measured emission 
line extinction for six volume limited luminosity bins. The 2-D bins extend 0.1 magnitudes 
in extinction and 2% in aperture fraction. The volume limited bins are defined as described 
in Figure [H The red contour shows the area where the 2-D bins contain at least 50 galaxies. 
The shading levels are described by the color bars on the left. 
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Fig. 13. — MC simulation results binned by Mr^o.i ior the 111,806 galaxies in the star forming 
class (diamonds) a r id th e 18,796 galaxies in the low S/N star forming class (triangles) of 
Brinchmann et al.l (120041 ). Each symbol represents a bin of 500 galaxies. The thin lines 
represent the upper and lower 95% confidence region measured for each star forming bin, 
while the thick lines do so for the low S/N bins. Top panel: The best fitting T values as a 
function of r-band luminosity. Lower panel: The values for each luminosity plotted on 
a log scale. Because the sample is dominated by the star forming class the results for these 
galaxies is very similar to the result for the whole sample in Figure [6l 
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Fig. 14. — Similar to Figure [S, this figure shows the value of as a function of the assumed 
/ ratio for the entire sample assuming F = 1.35. 
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Fig. 15. — Upper panel: The best fitting / ratio, assuming a universal Salpeter IMF, as a 
function of luminosity. The luminosity bins are volume limited and the same as in Figure [3 
Lower panel: The corresponding modified values as a function of luminosity. 
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Fig. 16. — Distribution in [g — r)o.i - log(Ha EW) space of the 329 galaxies in the volume 
limited M^,o.i = -17 bin (left) and the 8,049 galaxies in the M^,o.i = -23 bin (right). The 
contour levels and other descriptions are identical to Figure [TJ 
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Fig. 17. — Color and Ha EW values as a function of age for a model galaxy experiencing a 
burst (solid line). The model galaxy has solar metallicity, F = 1.35, and an exponentially 
decreasing SFH with r = 2.15 Gyr. The burst occurs at an age of 4.113 Gyr, lasts 250 Myr 
and has a strength of 10% of the total stellar mass. The age of the galaxy increases from 
the upper left to lower right before the bursts and the black dots appear on the track at 100 
Myr intervals. The dotted line shows the track for the model had a burst not occurred. 
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Fig. 18. — Best fitting IMF slope as a function of age for the model galaxy in Figure [T71 
which experiences a burst of star formation at an age of 4.113 Gyr. The dotted line indicates 
the underlying IMF model (F = 1.35). The F values have been boxcar smoothed by 20 Myr 
for clarity. 
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Fig. 19. — Distributions of best fitting parameters for the synthetic Mr,o.i = —17 data for 
models with Z = 0.01 and a single burst SFH. In all panels the dotted lines indicate the 95% 
errors determined from MC simulations. In the top two panels are results for T = 1.35. In 
the bottom two panels are results for F = 1.80. In the left column is the distribution of best 
fitting burst strengths given as a fraction of the total stellar mass formed. At right is the 
distribution of best fitting ages measured relative to the onset of the bursts. The bursts of 
star formation begin at and end at 200 Myr. For both IMFs 11 + 3/ — 5 of the 329 galaxies 
are best fit by models which have yet to experience their bursts and are not plotted. 
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Fig. 20. — Distributions of best fitting times measured from the beginning of a burst or gasp 
for the synthetic Mr,o.i = —17 data for models with Z = 0.01, F = 1.35 and multiple burst 
SFHs models described in Table [2l The values for each family of models is also shown. 
The dotted lines indicate the 95% confidence regions determined by MC simulations. 



-67- 




Fig. 21. — Histogram of the measured H5a values for the Mr,o.i = —17 (bold line) and 
Mr-,0.1 = ~23 (thin line) volume limited bins from Figure [16] expressed as a fraction of the 
total number of galaxies in each luminosity bin. The median uncertainty in H5a is 0.9 in both 
bins. The dotted lines are Gaussian profiles fitted to the distributions. For the M^^o.! = — 23 
bin (thin dots) the profile is centered on 3.4 with a = 1.8 and for M,.^o.i = bin (thick 
dots) the profile is centered on 6.1 with a = 0.9. 
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Fig. 22. — Model values of the H5a index as a function of age for four IMF models, F = 1.00 
(upper left), 1.40 (upper right), 1.70 (lower left), and 2.00 (lower right). Each dot represents 
a model value for single IMF, metallicity, SFH (as described in the models section) and age. 
The four tracks identifiable towards the left of each panel are each of differing metallicities 
with the lowest metallicity models having the largest US a values at young ages. The grey area 
indicates the range of B.6a spanned by the middle 90% of the volume limited M^^o.i = —23 
bin. The hashed area shows the middle 90% range for the volume limited M^ o.i = ~17 bin. 



